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CHAPTER  X 


GUIDE  TO  THE  FORTRAN  PROGRAM 


The  primary  functions  of  each  subroutine  as  well  as  the 
definitions  of  important  variables  are  given  in  the  following 
sections.  The  flags  and  conventions  which  are  unique  to  the 
HELP  code  are  also  listed  and  discussed. 

10.1  A  GENERAL  FLOW  DIAGRAM 

The  flow  diagram  given  in  Figure  10.1  indicates  the 
organization  of  the  code  and  describes  the  functions  of  the 
primary  subroutines  which  are  called  from  the  main 

driver  routine.  The  subroutines  which  in  turn  are  called 
by  each  of  the  primary  subroutines  are  also  indicated.  The 
calculational  cycle  has  six  stages:  (1)  the  calculation  of 
the  cell  pressures  and  the  time  step  (CDT) ;  (2)  the  periodic 

editing  of  cell  quantities  and  checking  for  energy  conservation 
(EDIT);  (3)  the  calculation  of  the  effects  of  deviator 
stresses  (SPHASE) ;  (4)  the  calculation  of  mass  transport 

terms  for  interface  cells  (INFACE);  and  (6)  the  calculation  of 
the  effects  of  transport  (TPHASE) .  For  calculations  which  ignore 
the  strength  effects,  the  third  stage  is  skipped.  For  calcula¬ 
tions  in  which  a  plug  is  generated,  another  stage  (PLGGEN)  ,  fol¬ 
lowing  SPHASE,  is  added  to  the  cycle.  The  rezone,  if  used,  is 
invoked  following  the  EDIT  stage. 


10.2  DESCRIPTION  OF  THE  SUBROUTINES 

The  primary  functions  of  each  HELP  subroutine  are  described 
below.  The  subroutines  are  listed  alphabetically. 

ADDENG 

ADDENG  is  called  only  if  one  of  the  material  packages  is 
a  high  explosive.  ADDENG  examines  the  detonation  time  associated 
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1)  Called  only  when  generating  problem  on  cycle 

2)  Called  only  when  detonating  a  high  explosive. 

3)  Called  only  whan  using  the  plugging  model. 

4)  Called  only  when  using  sliplines 


Figure  10. 1- -General  flow  diagram  of  the  HELP  code. 
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with  each  cell  containing  undetonated  high  explosives  to 
ascertain  if  it  is  time  for  it  to  be  detonated.  Once  it  is 
determined  that  cell  k  is  to  be  detonated,  DETIM(k)  is  set  to 
0,  which  allows  the  cell  to  become  active,  and  the  energy  of 
detonation  is  added  to  the  cell  and  system.  The  DETIM  array 
used  by  ADDENG  is  computed  in  DETIME  when  the  problem  is 
generated.  After  all  cells  containing  high  explosives  have 
been  detonated,  ADDENG  is  no  longer  called. 

ADDTCR 

Material  tracers  can  be  added  in  a  specified  region  at 
any  point  in  the  calculation  by  defining  input  flags  (NADD, 
NTRACR,  MINX,  MINY,  MAXX,  MAXY) .  When  NADD  >  0,  ADDTCR  is 
called  from  EDIT  every  cycle  which  is  a  multiple  of  NADD. 

ADDTCR  adds  new  tracers  between  existing  pairs  of  tracers  in 
the  region  defined  by  MINX,  MINY,  MAXX,  MAXY  but  does  not 
create  a  greater  density  of  tracers  than  that  specified  by 
NTRACR.  (See  Section  8.2.) 

CALFRC 

CALFRC  is  part  of  the  problem  generator  and  is  called  only 
in  the  initial  setup  of  a  problem.  Subroutine  CALFRC  follows 
each  tracer  string  through  the  grid  and,  for  each  cell  inter¬ 
sected  by  it,  determines  which  side  of  the  cell  the  string 
enters  (storing  this  information  in  the  variable  KODENT)  and 
through  which  side  the  string  leaves  (storing  this  information 
in  the  variable  KODLEV) .  The  X  and  Y  centimeter  coordinates 
of  the  entry  point  on  the  cell  boundary  are  stored  in  the 
first  location  of  the  arrays  XT  and  YT,  respectively.  The 
centimeter  coordinates  of  all  successive  tracer  particles 
which  lie  inside  the  cell  are  stored  in  succeeding 
locations  in  the  XT  and  YT  arrays.  When  CALFRC  determines 
the  string  has  left  the  cell,  the  centimeter  coordinates  of 
the  exit  point  are  determined  and  stored  in  the  next  available 
location  of  the  XT  and  YT  arrays.  There  can  be  no  more  than 
50  tracer  particles  which  define  the  position  of  the  material 


interface  in  the  mixed  cell.  CALFRC  defines  two  other 
variables,  NTST  and  NTND,  which  specify  the  initial  and 
final  material  tracer  particle  numbers,  respectively,  which 
lie  in  the  cell.  These  numbers  are  used  in  VOLFND  to  deter¬ 
mine  if  the  slipline  lies  in  the  cell. 

CARPS 

The  input  parameters  stored  in  the  Z-block  of  blank 
common  are  read  by  subroutine  CARDS.  The  format  of  these 
input  cards  is  described  in  Section  7.2.1. 

CDT 

A  principal  function  of  this  routine  is  to  compute  a 

time  step  which  ensures  stability  of  the  finite  difference 

equations.  This  is  done  by  finding  the  cell  with  the  minimum 

ration,  D/w.  Here  D  denotes  the  minimum  of  the  cell's 

dimensions  (cm)  and  u  denotes  the  sum  of  the  cell's 

maximum  velocity  component  and  its  sound  speed  (cm/sec) .  For 

an  ideal  gas  the  sound  speed  is  computed  as  /yP/p  and  for 

other  materials  by  the  approximate  relation  C  *  CQ  +  E”  Jp~ 

where  P  is  the  pressure  in  the  cell.  The  coefficient  ¥  is  an 

input  parameter  and  is  described  more  fully  in  Section  10.3 

(See  BBAR) .  CQ  is  defined  as  J  A/ pQ ,  where  A  is  the  bulk 

modulus  of  the  material.  The  values  of  C  for  nineteen 

o 

materials  are  given  in  a  DATA  statement  in  the  "included" 
element, COMDIM.  In  a  multimaterial  cell  the  sound  speed  is 
given  by  a  mass  weighted  average  of  the  sound  speeds  of  the 
materials  in  the  cell.  Each  cycle  CDT  prints  the  column  and 
row  (I,J)  of  the  cell  controlling  the  time  step  as  well  as  the 
maximum  sum  of  sound  speed  and  velocity  for  any  cell  in  the 
grid  (MAXCUV) ,  the  maximum  velocity  for  any  cell  in  the  grid 
(MAXUV) ,  and  the  velocity  and  pressure  cutoff  values  (UMIN,PMIN). 

Another  function  of  CDT  is  to  equilibrate  the  pressures 
of  materials  in  mixed  cells  by  using  an  iteration  scheme  that 
adjusts  the  material  densities.  A  detailed  discussion  of  this 
iteration  method  is  found  in  Section  4.6.1.  The  pressures  for 


pure  cells  are  also  updated  in  CDT  by  a  call  to  the  equation 
of  state  subroutine,  EQST .  The  density,  specific  internal 
energy  and  material  code  number  of  the  cell  (RHOW,  ENERGY,  N) 
are  passed  through  blank  common. 

When  a  high  explosive  is  being  detonated,  CDT  calls 
ADDENG  to  release  energy  in  cells  containing  H.E.  according 
to  the  detonation  time  (DETIM)  associated  with  each  cell. 

CMPRSN 

CMPRSN  computes  an  average  compression  (p/pQ)  of  the 
material(s)  in  a  cell.  These  compressions  are  used  as  weighting 
factors  by  HPHASE  in  the  definition  of  cell  boundary  pressures 
and  velocities  when  at  least  one  of  the  cells  contains  the 
free  surface  or  when  the  ratio  of  the  average  compressions  of 
the  two  cells  is  large  (>  CRATIO) .  (See  Section  2. 2. 2. 4.) 

COMDIM 

COMDIM  is  not  a  subroutine,  but  is  an  element  which  is 
inserted  in  place  of  the  statement  "INCLUDE  COMDIM"  in  most 
of  the  HELP  subroutines.  (This  terminology  is  unique  to 
UNIVAC  computers,  but  the  feature  is  offered  by  most  large 
computer  systems.) 

This  included  element  in  HELP  contains  dimension,  common 
and  equivalence  statements  and  two  data  statements  which  define 
material  constants,  pq  and  Cq,  for  various  solids  and  high 
explosives . 

COMPRS 

COMPRS  is  called  by  subroutine  REZONE  to  create  a  new  cell 
by  combining  two  adjacent  cells  into  one,  the  velocity  and 
energy  of  the  new  cell  being  defined  so  that  energy  and  momen¬ 
tum  are  conserved.  If  both  cells  are  pure,  then  the  new  cell 
will  also  be  pure.  If  one  or  both  contain  an  interface,  the 
new  cell  will  be  an  interface  cell.  (See  Section  8.1.) 


10-5 


BgJP 1  v  ,y  ryym  ^  » \  yj yu y^p  .■*  I  IJ1 MIJI  J.'J/l.M.1  »Mt»Pyf^PP^PWfP 


DETIME 

DETIME  defines  the  detonation  time  array,  DETIM,  used  in 
ADDENG.  DETIME  accomplishes  this  by  computing  the  time  re¬ 
quired  to  traverse  a  line  from  an  initiation  point  to  the 
center  of  a  cell,  plus  any  delay  time  associated  with  the 
initiation  point.  All  cells  in  the  path  must  contain  ex¬ 
plosives,  as  any  non-explosive  cell  in  the  path  is  treated 
as  an  obstacle  and  will  prevent  detonation  of  the  cell  by 
the  initiation  point.  The  X-  and  Y-coordinates ,  delay  time 
(used  by  primary  initiation  points  only),  and  material  package 
number  are  read  in  for  each  primary  and  secondary  initiation 
point.  (Secondary  initiation  points  are  defined  as  initiation 
sites  which  are  detonated  by  a  primary  initiation  point.) 

After  all  initiation  points  are  read  in,  the  X-  and  Y- 
coordinates  which  describe  the  probable  area  of  detonation 
(the  area  to  be  searched)  corresponding  to  each  of  the 
primary  and  secondary  initiation  points  are  read  in.  The 
detonation  time  is  computed  for  all  explosive  cells  within 
the  specified  area  to  which  there  is  an  unobstructed  path 
from  the  initiation  point.  See  Section  2.3.4 .1 .  The  input 
cards  are  described  in  Section  7.2.7  aim  in  Appendix  B. 

DMADJ 

DMADJ  adjusts  the  mass  transport  terms  of  interface  cells. 
On  the  first  subcycle  of  INFACE,  DMADJ  is  called  to  define  the 
mass  transport  terms  of  a  material  which  must  be  evacuated 
from  a  cell  because  its  interface  has  left  the  cell.  The  out- 
flux  terms  are  adjusted  (in  proportion  to  the  outflux  terms 
calculated  on  the  previous  cycle)  so  as  to  exactly  equal  the 
mass  remaining  in  the  cell.  The  influx  terms  are  set  to  zero. 

After  all  INFACE  subcycles  are  completed,  DMADJ  is  called 
again  to  check  that  the  computed  transport  terms  will  not  re¬ 
sult  in  a  negative  mass.  Also,  if  INFACE  is  being  subcycled 
(CYCMX  >  1),  the  evacuation  procedure  described  above  is 
repeated  since  an  interface  could  have  left  a  cell  on  a 
subsequent  subcycle.  (See  Section  4.4.) 


DMCALC 


DMCALC  computes  the  transport  terms  for  each  material 
across  the  four  boundaries  of  each  interface  cell.  DMCALC 
uses  the  fractional  cell  face  areas  (computed  by  subroutine 
FRACS)  associated  with  each  material  package  in  the  problem. 

In  this  way  the  position  of  the  material  interfaces  in  the  cell 
is  taken  into  account.  The  transport  terms  for  the  top,  right, 
bottom  and  left  boundaries  of  the  interface  cells  are  stored  by 
DMCALC  in  the  SAMPY,  SAMMP,  SAMMY,  and  SGAMC  arrays,  respect¬ 
ively.  TPHASE  performs  the  actual  transport  of  mass,  momentum 
and  energy  for  both  interface  and  pure  cells. 

EDIT 

The  periodic  printing  and  writing  on  the  restart  tape  are 
executed  by  subroutine  EDIT.  The  frequency  of  printing  and 
tape  dumps  is  controlled  by  input  parameters  described  in 
Section  7.2.1. 

On  every  print  cycle  EDIT  prints  the  mass,  total  energy, 
internal  energy,  kinetic  energy,  axial  and  radial  momentum, 
and  plastic  work  for  each  material  package  as  well  as  for  the 
entire  grid.  The  changes  in  energy  due  to  "evaporation"  and 
to  mass  lost  across  grid  boundaries  are  also  accounted  for. 

The  coordinates  of  the  material  tracers  circumscribing  each 
package  are  printed  in  cell  units.  Summary  graphs  of  the 
compression  (or  density),  pressure,  radial  velocity,  axial 
velocity  and  internal  energy  are  printed  if  the  input  param¬ 
eter,  MAPS,  is  non- zero.  And  finally,  the  pressure,  radial 
velocity,  axial  velocity,  internal  energy,  mass,  compression 
(or  density)  and  stress  deviators  for  all  cells  in  the 
active  grid  are  displayed  on  "long"  EDIT  prints  and  for  only 
the  cells  in  the  axis  column  on  "short"  EDIT  prints. 

EDIT  periodically  calls  ADDTCR  to  add  tracer  particles 
if  the  input  parameter  NADD  is  non-zero.  (See  Section  8.2.) 


Another  important  function  of  EDIT  is  to  check  the 
relative  error  between  the  sum  of  the  total  energy  of  all 
the  cells,  and  ETH,  the  theoretical  energy  of  the  grid.  If 
this  error  exceeds  a  limit  specified  by  the  input  variable 
DMIN,  EDIT  calls  ERROR  and  the  code  performs  an  error  exit. 

EDIT  also  determines  when  execution  should  stop  and 
sets  the  exit  flag,  WFLAGL. 

ENCHCK 

ENCHCK  is  called  by  the  main  routine,  HELP,  after  each 
phase  of  the  cycle  to  compute  the  relative  error  between  the 
total  grid  energy  and  ETH.  The  energy  of  the  cells  is  summed 
and  compared  to  the  theoretical  energy,  ETH,  which  is  the 
energy  at  time  T  =  0  plus  any  energy  added  to  the  grid 
(e.g.  when  an  explosive  is  being  detonated),  minus  any  energy 
subtracted  from  the  grid  (e.g.,  when  material  is  "evaporated" 
or  transported  out  of  the  grid) .  ENCHCK  prints  the  name  of 
the  phase  (TPHASE,  SPHASE,  HPHASE)  which  has  just  been 
completed,  along  with  ETH,  the  energy  sums  and  the  relative 
error.  Two  energy  sums  are  computed,  one  using  material 
velocities  (EMIX)  and  one  using  cell  velocities  (ESUM) .  The 
relative  error  is  based  on  the  EMIX  sum. 

These  prints  can  be  used  to  track  down  errors  which  re¬ 
sult  in  energy  not  being  conserved  by  noting  if  there  are 
consistent,  large  increases  in  the  relative  error  each  cycle 
following  the  execution  of  one  of  the  phases. 

ENDMV 

ENDMV  is  called  only  when  the  plugging  model  is  being 
used.  After  the  plug  is  completely  formed,  all  the  tracers 
on  the  vertical  edge  of  the  plug  are  moved  with  projectile 
and  plug  velocities  except  the  tracer  at  the  top,  which  is 
moved  with  target  velocities  only.  The  larger  axial 
velocities  in  the  plug  will  cause  the  tracers  below  the 
top  point  on  this  plug  surface  to . eventually  be  moved  above 
the  top  point  by  MOVTCR.  Subroutine  ENDMV  relocates 


these  tracers  at  the  position  of  the  top  point,  thereby 
accumulating  tracers  at  that  point  as  the  plug  moves  beyond 
the  back  surface  of  the  target. 

EOUT 

EOUT  keeps  track  of  the  changes  in  kinetic  and  internal 
energy  of  each  material  package  due  to  changes  in  the  theo¬ 
retical  grid  energy,  ETH.  EOUT  is  called  from  HPHASE,  TPHASE, 
and  SPHASE. 

EQST 

Given  a  material  code  number  as  well  as  a  density  and 
specific  internal  energy,  EQST  computes  a  pressure.  The 
equations  of  state  used  by  EQST  are  discussed  in  Section  2.3.1. 
The  material  constants  defined  in  DATA  statements  in  EQST  are 
listed  in  Tables  2.1  and  2.2.  CDT  calls  EQST  when  computing 
the  pressure  of  pure  cells  as  well  as  when  iterating  to 
equilibrate  material  pressures  in  multimaterial  cells. 

ERROR 

This  subroutine  is  called  in  case  certain  error  condi¬ 
tions  are  sensed  by  the  code.  ERROR  prints  a  message  identi¬ 
fying  the  general  location  within  the  subroutine  which  de¬ 
tected  the  error  condition  (see  error  message  #6  in  Section  9.1). 
It  lists  the  Z-block  variables  (the  first  150  words  of  blank 
common),  and  calls  EDIT  to  do  a  long'  print  and  a  tape  dump, 
and  then  stops  execution. 

FILGRD 

FILGRD  is  part  of  the  problem  generator  and  is  executed 
only  during  the  initial  setup  of  the  problem.  FILGRD  deter¬ 
mines  the  MFLAG  values  of  all  pure  cells  in  the  grid  except  for 
the  one  in  the  lower  left  corner,  which  the  user  defines. 

Given  the  initial  density,  velocity  components  and  internal 
energy  of  each  material  package,  FILGRD  defines  the  mass, 
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velocity  components  and  specific  internal  energy  of  each  cell 
in  the  grid.  The  mass  of  a  pure  cell  is  based  on  the  cell's 
MFLAG  value  which  indicates  which  package  the  cell  belongs  to. 
If  the  cell  is  an  interface  cell,  the  mass  of  each  material 
in  the  cell  is  computed  by  FILGRD  using  the  partial  volumes 
computed  by  VOLFND. 

FLGSET 

FLGSET  is  called  on  each  subcycle  of  INFACE,  after  the 
fractional  cell  face  area  terms  are  computed.  Based  on  the 
updated  area  terms,  FLGSET  determines  which  (if  any)  inter¬ 
face  cells  have  become  pure  (lost  all  interfaces),  and  which 
(if  any)  have  lost  some,  but  not  all,  interfaces.  In  either 
case  one  or  more  materials  must  be  evacuated  from  the  cell. 
FLGSET  implements  the  following  conventions  to  denote  the 
occurrence  of  these  events:  (1)  the  density  RHO(n,m)  of  any 
material (s)  to  be  evacuated  is  set  to  zero;  (2)  the  MFLAG  of 
any  interface  cell  that  has  become  pure  is  made  negative. 
However,  any  cell  that  has  become  pure  is  still  handled  as 
an  interface  cell  until  the  end  of  the  cycle,  when  its  MFLAG 
is  changed  in  TPHASE. 

FLGSET  also  checks  that  the  sum  of  the  fractional  cell 
face  areas  for  each  interface  cell  boundary  equals  the  total 
area  of  that  boundary.  If  the  sum  differs  from  the  total 
area  by  more  than  1%  an  error  message  is  printed  (see  error 
message  #30  in  Section  9.1).  The  error  is  usually  the  result 
of  interfaces  crossing  because  the  tracer  particles  have  be¬ 
come  too  sparse. 

FRACS 

The  straight  line  segments  which  connect  the  tracer  par¬ 
ticles  that  circumscribe  each  material  package  are  used  by 
FRACS  to  determine  where  the  interfaces  cut  cell  boundaries 
and  to  compute  the  fractional  cell  face  areas.  Section  4.3 
gives  a  detailed  description  of  the  logic  and  conventions 
employed  by  FRACS. 
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HELP 

The  overall  cycling  of  the  calculation  is  controlled  by 
HELP,  as  shown  in  Figure  10.1.  HELP  is  the  main  routine  of 
the  code  and  calls  INPUT,  CDT,  EDIT,  REZONE  (if  it  is  a  rezone 
cycle),  SPHASE,  PLGGEN  (if  the  plugging  option  is  used), 

HPHASE,  INFACE  and  TPHASE  in  that  order.  When  the  user  desires 
to  see  the  effects  of  each  phase  of  the  calculational  cycle, 

HELP  will  respond  to  the  input  parameter  INTER  and,  on  print 
cycles,  call  EDIT  after  SPHASE  and  HPHASE  as  well  as  after  CDT 
(see  Section  9.2).  HELP  also  calls  EXIT  on  the  normal 
cessation  of  the  calculation  when  WFLAGL  >  0. 

HPHASE 

The  effect  of  the  pressure  gradients  in  updating  the 
velocities  and  the  internal  energies  is  computed  in  HPHASE. 

(The  numerical  method  is  described  in  detail  in  Section  2.2.2.) 

INFACE 

INFACE  controls  the  subcycling  of  FRACS,  FLGSET,  DMADJ, 

DMCALC,  and  MOVTCR.  the  subroutines  which  compute  the  mass  trans¬ 
port  terms  of  interface  cells  and  which  advance  the  tracer 
particles.  If  sliplines  are  used  (NOSLIP  =  0)  INFACE  calls 
THETAS,  after  the  subcycles  are  completed,  to  compute  the 
angle  of  the  slipline  across  each  slipline  cell. 

INPUT 

Instructions  for  running  problems  are  interpreted  by  sub¬ 
routine  INPUT,  which  can  either  start  or  restart  a  calculation. 
It  calls  CARDS  to  read  the  first  part  of  the  input  deck,  and, 
on  cycle  0,  calls  SETUP  to  generate  the  initial  conditions  of 
the  calculation.  When  generating  a  calculation,  INPUT  rewinds 
and  reads  the  restart  file  after  SETUP  writes  the  cycle  0  dump, 
and,  when  restarting  a  calculation,  INPUT  finds  and  reads  the 
restart  dump. 

MAP 

This  subroutine  is  called  by  EDIT  when  the  input  parameter 
MAPS  is  non-zero.  It  displays  the  properties  of  each  cell  in 
the  active  grid  using  an  alphabetic  scale.  The  user  obtains 
symbolic  maps  of  the  compression,  pressure,  radial  and  axial 
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velocities,  and  specific  internal  energy  by  setting  MAPS  =  1. 

To  get  a  density  map  instead  of  a  compression  map,  the  user 
sets  MAPS  =  2.  Asterisks  are  displayed  for  the  interface  cells 
in  the  density  and  compression  maps.  The  scale  of  each  map  is 
adjusted  according  to  the  current  minimum  and  maximum  values 
of  each  quantity. 

MOVTCR 

The  material  and  passive  tracer  particles  are  moved  with 
a  local  velocity  field  in  MOVTCR.  A  description  of  the  method 
for  moving  tracers  is  given  in  Section  4.1.2. 

NEWFLG 

NEWFLG  determines  the  MFLAG  value  of  new  interface  cells 
(see  Section  4.2.1). 

NEWMIX 

NEWMIX  is  called  from  FRACS  to  define  variables  for  a  new 
interface  cell.  If  a  cell  becomes  an  interface  cell  after  the 
first  subcycle  of  INFACE,  NEWMIX  calls  DMCALC  to  make  the  def¬ 
inition  of  the  mass  transport  terms  for  that  cell  consistent 
with  the  number  of  subcycles  already  completed.  (See  Sections 
4.2.2  and  4.2.3.) 

NEWRHO 

When  the  material  n  interface  first  enters  a  cell,  INFACE 
calls  subroutine  NEWRHO  to  define  the  density  of  material  n  for 
that  cell  by  looking  at  the  density  of  material  n  in  the  neigh¬ 
boring  cells  closest  to  the  interface.  (See  Section  4.2.2.) 


PLGADD 

PLGADD  is  called  by  SETUP  only  when  the  plugging  model  is 
being  used.  It  adds  target  and  plug  material  tracer  particles 
at  the  corner  where  the  plug  surface  will  begin.  Those  addi¬ 
tional  points  will  be  used  to  define  the  vertical  edge  of  the 
plug  as  it  evolves.  This  edge  will  be  a  slipline  and  the  plug 
will  be  package  2  and  a  "slave",  so  the  number  of  tracers 
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added  at  the  corner  is  defined  by  the  plug  slipline  end¬ 
points,  NBGSD(2)  and  NF.HDSD(2),  which  in  turn  are  defined  by  the 
user  in  the  INPUT  deck.  (See  Section  6.1.) 

PLGALF 

PLGALF  computes  the  angle  of  maximum  shearing  stress  for 
cells  in  the  plugging  region  of  the  target,  given  their  deviator 
strain  rates.  PLGALF  is  called  from  SPHASE  only  if  the  plugging 
model  is  being  used. 

PLGGEN 

At  the  end  of  SPHASE,  when  the  plugging  model  is  activated, 
PLGGEN  is  called  to  update  the  specific  plastic  work  of  the 
material  associated  with  each  passive  tracer  in  the  plugging 
region  of  the  target.  If  the  criterion  for  extending  the  plug 
surface  is  satisfied,  PLGTCR,  PLGVOL  and  PLGMAS  are  called 
to  enlarge  the  plug  package.  (See  Section  6.2.1.) 

PLGMAS 

PLGMAS  converts  target  material  (package  3)  into  plug 
material  (package  2)  by  redefining  the  RHO,  SIE,  and  XMASS 
arrays.  PLGMAS  also  updates  the  plug  and  slipline  endpoints 
(SLPNDX ,  SLPNDY) .  (See  Section  6.2.4.) 

PLGTCR 

When  the  slipline  is  advanced,  PLGTCR  moves  the  tracers 
at  the  top  of  the  plug  package  into  the  row  containing  the 
new  end  point  of  the  slipline.  When  the  last  row  of  the 
target  fails,  the  tracers  at  the  top  of  the  plug  are  made 
coincident  with  the  target  tracers  out  to  the  point  at 
which  the  slipline  intercepts  the  top  of  the  target. 

Two  plug  tracers  are  left  at  that  intercept  point.  One  moves 
out  with  the  plug,  the  other  remains  attached  to  the  target 
and  is  the  end  point  of  the  slipline.  The  target  tracers  are 
also  redefined.  The  target  package  is  no  longer  attached  to 
the  axis,  but  begins  at  the  top  slipline  endpoint.  (The  tracers 


between  the  axis  and  the  slipline  are  plug  and  free  surface 
tracers.)  The  indices  of  the  tracers  of  each  package  that 
define  the  endpoints  of  the  slipline (NBGSD^ ,  NBGMD^,  NENDSD^, 
NENDMD^,  i  =  i,  NMAT)  are  changed  as  necessary  at  this  time. 

(See  Section  6.2.2.) 

PLGVOL 

When  the  plug  is  extended  through  another  row  of  cells, 

PLGVOL  is  called  to  compute  the  partial  cell  volume  of  the 
plug  material  in  the  cell  containing  the  extended  plug 
surface.  (See  Section  6. 2.4. 2.) 

PLUGUV 

Once  the  plug  has  begun  to  evolve,  the  motion  of  the 
tracers  on  the  plug  interface  is  determined  by  special  pre¬ 
scriptions  in  PLUGUV.  These  prescriptions  are  described  in 
Section  6.2.3. 

PTSAV 

The  intercepts  between  the  slipline  and  grid  lines  are 
stored  by  subroutine  PTSAV  and  used  by  subroutine  THETAS  to 
compute  the  angle  of  the  slipline  across  each  slipline  cell. 

PTSAV  is  called  by  FRACS  when  a  pair  of  master  package  tracers 
along  the  slipline  are  being  used  to  compute  a  fractional  cell 
face  area.  FRACS  computes  the  intercept  (in  cm.)  and  PTSAV 
stores  it  in  special  arrays,  XINT  and  YINT. 

PTSAV  uses  a  pointer  array,  MSLD,  to  associate  the  MFLAG 
and  the  i  and  j  values  of  a  tell  with  the  intercept  arrays, 

XINT  and  YINT,  as  follows: 

MSLD  (ms)  =  MFLAG  (k)*1000000  +  i*1000  +  j 

XINT  (l,ms)  -*•  x-coordinate  of  the  intercept  at  top  of  cell  (k,i,j) 
XINT  (2, ms)  -*  x-coordinate  of  the  intercept  at  bottom  of  cell  (k,i, 
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YINT  (l,ms)  -*•  y-coordinate  of  the  intercept  at  right  of  cell  (k,i,j) 

YINT  (2, ms)  -*•  y-coordinate  of  the  intercept  at  left  of  cell  (k,i,j). 

The  XINT,  YINT  arrays  are  initialized  on  each  subcycle  of 
INFACE. 

PTSAV  also  sets  THETA  (m)  =  100  to  indicate  that  the  mixed 
cell  k  is  cut  by  the  slipline,  where  m  =  MFLAG  (k)  -  100. 

REZONE 

To  add  more  material  and  diminish  the  number  of  zones  in 
the  active  grid  the  REZONE  routine  combines  either  two  cells 
(if  rezoning  in  only  one  direction)  or  four  cells  (if  rezoning 
in  both  directions)  into  one,  and  adds  Columns  and/or  rows  to 
maintain  an  IMAX  by  JMAX  grid.  (See  Section  8.1.) 

REZONE  is  called  from  HELP  only  when  certain  input  param¬ 
eters  (NUMREZ,  REZ,  IEXTX,  JEXTY)  are  set.  It  can  be  "forced"  on 
the  first  cycle  of  a  restart  run  by  setting  REZ=1  in  the  input 
deck,  or  it  can  be  "triggered"  by  TPHASE  which  sets  REZ-1  when 
a  signal  reaches  the  edge  of  the  grid. 

RNDOFF 

RNDOFF  is  called  from  DMADJ  to  insure  that  the  sum  of  the 
mass  transport  terms  of  a  given  material  exactly  equals  the 
mass  of  that  material  remaining  in  the  cell.  RNDOFF  is  called 
when  a  material  is  being  evacuated  or  its  transport  terms  are 
being  adjusted  to  prevent  a  negative  mass.  (See  Section  4.4.) 

SETUP 

SETUP  defines  several  problem  constants  and  reads  the 
cards  defining  the  DX  and  DY  arrays,  the  initial  conditions 
and  strength  properties  of  the  material  packages,  and  the 
parameters  defining  the  slipline.  SETUP  also  generates  the 
passive  tracer  particles  and  writes  the  cycle  0  tape  dump. 
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SETUPA 

SETUPA  defines  the  X,  Y  and  TAU  arrays  and  converts  the 
tracer  particle  positions  from  centimeter  to  cell  units.  It 
moves  tracer  particles  that  are  close  to  a  grid  bound¬ 
ary  onto  that  boundary,  and  forces  positions  of  corresponding 
tracers  from  different  packages  to  be  exactly  equal.  SETUPA 
also  initializes  some  of  the  mixed  cell  arrays  and  calculates 
the  initial  value  of  ETH,  the  total  theoretical  energ)  of  the 
material  in  the  grid.  CALFRC  and  FILGRD  are  called  from 
SETUPA  to  define  the  properties  of  each  cell  in  the  grid.  For 
plugging  calculations  PLGADD  is  called  from  SETUPA  to  generate 
additional  tracer  particles 

SPHASE 

In  subroutine  SPHASE,  the  deviator  stresses  acting  on  each 
cell  face  and  the  hoop  stress  are  determined,  and  the  result¬ 
ing  velocity  and  internal  energy  increments  are  computed. 
(Details  are  given  in  Sections  2.2.1  and  2.3.2,)  SPHASE  calls 
STRNG  to  compute  the  strength  of  the  material  in  a  cell.  The 
effects  of  strength  can  be  omitted  and  SPHASE  bypassed  by 
setting  CYCPH3  «=  -  1. 

STRNG 

The  yield  strength  of  the  material  in  a  cell  is  computed 
by  STRNG.  (See  Section  2.3.2.)  The  strength  of  a  multi¬ 
material  cell  is  a  volume  weighted  average  of  the  strength  of 
its  constituents.  The  strength  constants  for  each  material 
are  defined  by  input  cards  when  the  problem  is  generated. 

STRNG  is  called  from  SPHASE. 

THETAS 

For  calculations  which  use  the  slipline  capability  of  the 
code  (NOSLIP=0),  THETAS  is  called  to  compute  the  angle  of  the 
slipline  in  each  slipline  cell  for  which  PTSAV  set  THETA (m) =100 . 
THETAS  uses  the  intercepts  saved  by  PTSAV  to  determine  the 
angle  of  the  slipline  across  a  given  slip  cell.  THETAS  uses 
the  first  two  non-zero  cell  boundary  intercepts  in  this  order: 


top,  bottom,  right,  and  left.  Therefore,  if  the  slipline 
cuts  the  top,  left, right,  and  bottom  sides  of  a  cell,  THETAS 
uses  only  the  top  and  bottom  intercepts  to  define  an  average 
angle  of  the  slipline  through  that  cell. 

TPHASE 

Mass  transport  and  the  associated  transport  of  momentum 
and  energy  are  accounted  for  in  TPHASE.  (The  numerical 
method  is  described  in  detail  in  Section  2.2.3.)  Before 
printing  a  symbolic  map  which  displays  the  material  package 
numbers  and  interface  cell  locations,  TPHASE  redefines  the 
flags  of  cells  that  have  become  pure. 

TSETUP 

TSETUP  reads  cards  defining  tracer  particle  locations. 
TSETUP  generates,  in  centimeter  coordinates,  tracer  particle 
positions  which  are  on  straight  line  segments  or  on  arcs  of 
circles  or  ellipses.  (See  Section  7.2.5.) 

UC 

Given  the  magnitude  of  the  velocity  components  normal  and 
tangential  to  the  slipline  (WN,  WT) ,  subroutine  UC  computes  the 
radial  component  of  velocity  in  the  grid  coordinate  system. 

UVCALC 

UVCALC  is  called  by  HPHASE  to  update  the  velocities  of 
materials  in  slipline  cells.  The  new  material  velocities  are 
determined  by  solving  the  four  simultaneous  linear  equations 
described  in  Section  5.2.  If  a  cell  does  not  contain  a 
slipline  (THETA(m)  <  0),  HPHASE  sets  the  material  velocities 
(US,  VS)  equal  to  the  updated  cell-centered  velocities  (U,  V). 

UVMOD 

After  material  and  cell  velocities  are  modified  in 
TPHASE  due  to  momentum  transport,  UVMOD  is  called  to  insure 
that  the  material  velocities  of  slipline  cells  satisfy  the 
requirement  that  the  master  and  slave  velocity  components 


normal  to  the  slipline  are  equal  and  that  the  total  momentum 
of  the  cell  is  conserved.  (See  Section  5.4.)  UVMOD  is  not 
called  if  sliplines  are  not  used  (NOSLIP  >  0). 

VC 

Given  the  magnitude  of  the  velocity  components  normal 
and  tangential  to  the  slipline  (WN,WT),  subroutine  VC  computes 
the  axial  component  of  velocity  in  the  grid  coordinate  system. 

VDCLOS 

Given  the  second  index  (NVRTEX)  of  the  free  surface  tracer 
particle  which  is  the  vertex  of  a  void  closing  region,  VDCLOS 
applies  three  criteria  to  determine  if  the  void  should  be  closed 
up  to  the  next  free  surface  tracer.  (See  Section  8.4.) 

If  one  of  the  criteria  is  met,  VDCLOS  moves  two  material 
tracers  and  one  free  surface  tracer  to  a  specified  point  on 
the  free  surface.  The  specified  point  then  becomes  the  new 
vertex  of  the  void  closing  region.  Two  free  surface  tracers 
are  removed,  and,  in  the  region  where  the  void  was  closed,  the 
interface  becomes  only  a  material  interface  and  is  no  longer 
a  free  surface. 

VOLFND 

VOLFND  is  part  of  the  problem  generator  and  is  called 
only  in  the  initial  setup  of  a  problem.  VOLFND  defines  the 
volume  and  surface  areas  of  a  given  material  package  in  a 
mixed  cell.  VOLFND  accomplishes  this  by  using  the  XT,  YT 
arrays  which  define  the  tracer  particles  belonging  to  that 
package  and  located  in  that  cell.  VOLFND  adds  tracer 
particles  at  appropriate  corners  of  the  cell,  as  well  as 
one  having  the  same  coordinates  as  the  first  tracer  in  the 
XT,  YT  array.  The  arrays,  XT,  YT,  now  define  an  enclosed 
volume  in  such  a  way  that,  as  one  proceeds  between  two  consec¬ 
utive  tracers,  the  given  material  package  is  to  the  left. 

VOLFND  calculates  the  partial  volume  and  fractional  surface 
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areas  for  that  material  through  the  use  of  the  theorems  of 
Pappus.  For  the  volume  calculation,  the  theorem  states  the 
volume  generated  by  a  plane  area  that  is  rotated  about  a 
line  that  lies  in  its  plane  but  does  not  intersect  the  area, 
is  equal  to  the  product  of  the  area  and  the  distance  traveled 
by  its  center  of  gravity.  For  the  surface  area  calculation, 
the  theorem  states  the  surface  area  generated  by  a  line  that 
is  rotated  about  a  line  that  lies  in  its  plane  but  does  not 
intersect  the  line,  is  equal  to  the  product  of  the  length  of 
the  line  and  the  distance  traveled  by  its  center  of  gravity. 

The  partial  volume  is  stored  in  the  VOLM  array.  The  fract¬ 
ional  surface  areas  are  stored  in  FRACRT  and  FRACTP. 

VOLFND  also  determines  if  the  slipline  lies  in  the  mixed 
cell  and,  if  so,  calculates  the  angle  THETA  the  slipline 
makes  with  the  X-axis. 

WN 

Given  radial  and  axial  velocity  components  and  the  angle 
of  the  slipline,  WN  computes  the  magnitude  of  the  velocity 
component  normal  to  the  slipline. 

WT 

Given  radial  and  axial  velocity  components  and  the  angle 
of  the  slipline,  WT  computes  the  magnitude  of  the  velocity 
component  tangent  to  the  slipline. 

XCTOP 

XCTOP  converts  the  x-coordinate  of  a  point  from  cell  to 
centimeter  units. 

XPTOC 

XPTOC  converts  the  x-coordinate  of  a  point  from  centimeter 
to  cell  units. 

YCTOP 

YCTOP  converts  the  y-coordinate  of  a  point  from  cell  to 
centimeter  units. 


YPTOC 


YPTOC  converts  the  y-coordinate  of  a  point  from  centimeter 
to  cell  units. 

10.3  DICTIONARY  OF  IMPORTANT  VARIABLES 

This  section  includes  a  description  of  the  use,  units, 
dimension  and  location  of  all  variables  in  common  blocks,  as 
well  as  many  variables  local  to  the  subroutines.  The  following 
conventions  are  used  in  the  dictionary  in  describing  the 
storage  location  of  the  variables: 

(NAME)  The  variable  is  local  to  subroutine  NAME. 

NAME  The  variable  is  in  common  block  NAME. 

B. C.  The  variable  is  in  Blank  Common  or_  equivalenced 

to  a  variable  in  Blank  Common. 

=  Z(N)  The  variable  is  equivalenced  to  a  member  of  the 

Z-array,  the  first  array  in  Blank  Common.  Most  of 
these  variables  are  used  in  generating  and  restart¬ 
ing  problems. 

---  The  variable  is  used  as  a  calling  argument. 

In  describing  the  dimensions  of  the  variables  the  conventions 
are : 

---  The  variable  is  not  dimensioned. 

(30)  The  array  is  always  dimensioned  30. 

(NVAR)  The  dimension  of  the  array  must  be  greater  than  or 

equal  to  the  value  of  NVAR,  where  NVAR  is  determined 
by  the  user  when  the  calculation  is  generated.  (See 
Appendix  A  for  procedures  for  dimensioning  the  HELP 
arrays . ) 


Variable 

Name 

Location 

Dimens  ion 
of  Array 

Units 

Definition 

A 

(EQST) 

_ . . . 

— 

"a"  in  the  thermal 

pressure  term  of 
the  Tillotson 


equation  of  state. 
(See  Section  2. 3. 1.2) 
Also,  "a"  in  the 
JWL  equation  of 
state  for  detonated 
explosives.  (See 
Section  2. 3. 1.1) 


AIX 

B.C. 

(KMAX) 

ergs/g 

Specific  internal 
energy  in  a  cell. 

ALE 

! 

(MAP) 

(41) 

This  array  has 
alphabetic  char¬ 
acters  for  posi¬ 
tive  valuer  in 
the  density  (or 
compression) , 
velocity*  pressure 
and  energy  maps. 
Defined  in  a  DATA 
statement . 

ALFA 

(EQST) 

— 

— 

"aM  in  the 

Tillotson  equation 
of  state  for  hot, 
expanded  materials 
(see  Section 
2 . 3 . 1 . 2).  "a"  in  the 
JWL  equation  of 
state  for  a 
detonated  high 
explosive.  (See 
Section  2. 3. 1.1) 

ALPHA  =Z (7 5)  ---  radians  In  PLGALF,  ALPHA 

is  the  current 
angle  of  maximum 
shearing  stress 
of  a  given  cell. 
Used  only  in 
plugging  calcu¬ 
lations.  In 
PLGTCR,  PLGVOL 
and  PLGMAS,  it  is 
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Variable 

Name 


Location 


Dimension 
of  Array 


Units 


ALPHSV 


(PLGGEN) 


radians 


ALPSV 


PLSTC 


(IPLGRT, 

IPLGTP- 

IPLGBT+1) 


radians 


AMDM 


B.C. 


(NMAT) 


AMMP 


(TPUASE) 


AMMU 


(TPHASE) 


g-cm/sec 
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Definition 


the  angle  at 
which  the  plug 
edge  is  being 
extended . 


The  angle  of  maxi¬ 
mum  shearing  stress 
of  material  (assoc¬ 
iated  with  a  pas¬ 
sive  tracer)  which 
has  just  satisfied 
the  plastic  work 
criterion  and  is 
near  the  tip  of 
the  plug. 


The  angle  of 
maximum  shearing 
stress  of  each 
cell  in  the  plug¬ 
ging  region  of 
the  target  is 
saved  in  the 
ALPSV  array. 


The  tensile 
failure  criterion 
is  based  on 
material  expansion. 
If  p/p  <  AMDM,  the 
materi  Si  cannot 
support  hydrostatic 
tensions  or  devi¬ 
ator  stresses.  (See 
Section  2.3.3) 


Mass  transported 
across  the  right 
boundary  of  a  cell. 
(See  Section 
10.4.2) 


Radial  momentum 
transported  across 


the  bottom  bound¬ 
ary  of  a  cell. 

(See  Section 
10.4.2) 

AMMV  (TPHASE)  -  g-cm/sec  Axial  momentum 

transported 
across  the  bottom 
boundary  of  a 
cell.  (See 
Section  10.4.2) 

AMMY  (TPHASE)  -  g  Mass  transported 

across  the 
bottom  boundary 
of  a  cell.  (See 
Section  10.4.2) 

AMPY  (TPHASE)  -  g  Mass  transported 

across  the  top  bound 

t  ary  of  a  cell. 

(See  Section 
10.4.2) 

AMUR  (TPHASE)  -  g-cm/sec  Radial  momentum 

transported 
across  the  right 
boundary  of  a 
cell.  (See 
Section  10.4.2) 

AMUT  (TPHASE)  -  g-cm/sec-.  Radial  momentum 

transported  across 
the  top  boundary  of 
a  cell.  (See 
Section  10.4.2) 

AMVR  (TPHASE)  -  g-cm/sec  Axial  momentum 

transported  across 
the  right  boundary 
of  a  cell.  (See 
Section  10.4.2) 

AMVT  (TPHASE)  -  g-cm/sec  Axial  momentum 

transported  across 
the  top  boundary  of 
a  cell.  (See 
Section  10.4.2) 


Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

AMX 

B.C. 

( KMAX ) 

g 

Total  mass  in  a 
cell . 

AN 

(PLGGEN) 

The  number  of 
passive  tracers 
near  the  tip  of 
the  plug  whose 
material  has 
passed  the  plas¬ 
tic  work  criterion 
on  a  given  cycle. 

B 

(EQST) 

"b"  in  the  ther¬ 
mal  pressure  term 
of  the  Tillotson 
equation  of  state. 
(See  Section 

2. 3. 1.2) 

BBAR 

*  Z(149) 

— 

— 

"E"  used  in  CDT.  An 

INPUT  parameter 
used  in  local 
sound-speed  calcu¬ 
lation  for  all 
materials  other 
than  an  ideal 
gas.  (Local 
sound-speed  for 
material  i  in  cell 
K  is  approximated 
as  (C0i)*ff’/PTKT. 
where  the 
coefficient  B  is 
obtained  by  de¬ 
termining  a 
typical  slope 
for  the  isentropes 
in  Ref. 11  and 
using  the  relation 
C  =  VZ-TP/dV  to 
evaluate  B  at  a 
particular  point. 

BBOUND  =Z(74)  — -  ergs  Calculated  in 

SPHASE.  Printed 
in  EDIT  under 
"Elastic  Plastic 
Work."  Total 
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Variable 

Name 


Location 


Dimension 
of  Array 


Units 


BETA 


(EQST) 


BOTM 


'Z  (39) 


BOTMU 


-Z(64) 


g-cm/sec 


BOTMV 


55  Z  ( 4  0  ) 


g-cm/ sec 


CAPA 


(EQST) 


dynes/cm' 


CAPB 


(EQST) 


dynes/cm 1 
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Definition 


work  done  by  the 
elastic  and 
plastic  stresses, 


"p"  in  the  Tillot- 
son  equation  of 
state  for  hot, 
expanded  materials. 
(See  Section 
2. 3. 1.2).  "3"  in 
the  JWL  equation. 
(See  Section 
2. 3. 1.1) 


Calculated  in 
T PHASE.  Printed 
in  EDIT.  Total 
mass  transported 
across  bottom  of 
grid. 


Calculated  in 
TPHASE.  Printed 
in  EDIT.  Total 
radial  momentum 
transported  across 
bottom  of  grid. 


Calculated  in 
TPHASE.  Printed 
in  EDIT.  Total 
axial  momentum 
transported  across 
bottom  of  grid. 


"A"  in  the  mechan¬ 
ical  pressure  term 
of  the  Tillotson 
equation  of  state 
and  of  the  JWL 
equation  of  state 
for  detonated 
explosives.  (See 
Sections  2. 3. 1.1 
and  2 . 3 . 1 . 2) 


"B"  in  the  mechan¬ 
ical  pressure  term 
of  the  Tillotson 
and  JWL  equations 
of  state.  (See 
Sections  2.3. 1.1 
and  2 . 3 . 1 . 2) 


eS_ 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


CKE 


CNAUT 


COMPC 


CRATIO 


CSQR 


(HPHASE) 


MXCELL 


(CMPRSN) 


=  Z (148) 


B.C. 


(30) 


(NMAT) 


ergs  The  kinetic  energy 

of  a  slipline  cell 
after  the  material 
velocities  are  up¬ 
dated  by  UVCALC . 

cm/sec  "C  "  used  in  CDT . 

Approximate  sound- 
speed  of  nineteen 
materials  defined 
in  a  DATA  State¬ 
ment  in  COMDIM,  the 
"included"  element . 

f  ESCAPA. 

r  .  -  /  _ I 

Loi  /  RHOZ. 

1 

*  >Vp0i 

-  The  compression  of 

the  material  in 
cell  k.  Used  in 
HPHASE  to  weight 
pressures  and 
velocities,  when 
necessary,  in 
computing  cell 
boundary  values. (See 
Section  2. 2. 2. 4) 

-  The  cell  boundary 

pressures  and 
velocities  are 
compression  weighted 
in  HPHASE  if  the 
ratio  of  the  com¬ 
pressions  of  the 
two  cells  is 
greater  than 
CRATIO,  an  input 
parameter.  (See 
Section  2. 2. 2. 4) 

dynes-cm/g  The  constant 

energy  compressi¬ 
bility  of  a 
material,  computed 
and  used  in  CDT 
in  the  pressure 
iteration.  (See 
Section  4. 6. 1.1) 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


CVIS  =Z  (2 7)  -  -  INPUT  parameter. 

Used  to  describe 
the  bottom  bound¬ 
ary  condition.  The 
bottom  grid  boundary 
is  transmittive 
when  CVIS  =  -1, 
reflective  when 
CVIS  =  0. 

CYC  B.C.  -  -  A  flag  used  in 

DMCALC,  DMADJ  and 
set  in  INFACE,  FRACS, 
SETUPA  and  NEWMIX. 

In  DMCALC  it  is 
used  as  a  factor 
in  the  computa¬ 
tion  of  the  mass 
transport  terms 
and  is  greater 
than  1  only 
when  NEWMIX  calls 
DMCALC  to  make  up 
subcycles  for 
cells  which  be¬ 
come  interface 
cells  after  the 
first  subcycle  of 
INFACE.  In  DMADJ, 
CYC  =  0  means  the 
mass  transport 
terms  have  not 
been  updated,  and 
CYC  =  1  means 
all  subcycles  of 
INFACE  have  been 
completed  and  the 
mass  transport 
terms  are  updated. 
(See  Section  4.4) 

CYCLE  »Z(2)  -  -  Used  in  INPUT,  MAP, 

SETUP,  CDT ,  EDIT. 
Cycle  number  (an 
integer  value  in 
floating  point 
form) . 


Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

CYCMX 

=  Z (69) 

INPUT  parameter. 
Equals  number  of 
passes  through 

INFACE  on  each 
cycle.  Suggested 
minimum  of  2  and 
maximum  of  8. 

Used  to  minimize 
transport  noise 
near  interfaces. 

CYCPH3 

aZ (7Q) 

INPUT  parameter. 
Number  of  times 
to  subcycle  SPHASE. 
If  it  is  set  to  -1, 
SPHASE  is  omitted. 

CZERO 

B.C. 

(NMAT) 

dynes/cm2 

Value  of  Y  for  cal¬ 
culation  of  material 
yield  strength. 
Defined  by  input 
cards  for  each 
material  package 
in  the  grid.  Used 
in  STRNG .  (See 
Section  2. 3. 2. 3) 

DELEB 

(TPHASE) 

ergs 

The  total  energy 
transported  across 
the  bottom 
boundary  of  a 
cell.  (See 

Section  10.4.2) 

DELER 

(TPHASE) 

ergs 

The  total  energy 
transported  across 
the  right 
boundary  of  a 
cell.  (See 

Section  10.4.2) 

DELET 

(TPHASE ) 

ergs 

The  total  energy 
transported  across 
the  top  boundary 
of  a  cell.  (See 
Section  10.4.2) 

DELM 

(TPHASE) 

.  _ . . 

g 

Total  change  in 
mass  of  a  cell. 

I 


Variable 

Name 

Location 

Dimens  ion 
of  Array 

Units 

Definition 

DENGY 

(ADDENG) 

(30) 

ergs/g 

The  specific  interna 
energy  released  by 
a  high  explosive, 
defined  in  a  DATA 
statement  in  ADDENG 
for  four  explosives. 

DETIM 

B.C. 

(KMAX) 

sec 

The  time  for  the  HE 
detonation  front 
to  reach  the  center 

of  a  cell.  If  cell 
K  does  not  con¬ 
tain  high  explo¬ 
sive,  or  if  the 
explosive  has  been 
detonated,  DETIM(K) 
=0.  If  the 
problem  does  not 
involve  an  explo¬ 
sive  the  DETIM 
array  can  be 
dimensioned  1. 

DISTX  (MOVTCR)  -  cm  The  distance  in 

the  x-direction  a 
tracer  moves  on  a 
given  subcycle  of 
INFACE. 

DISTY  (MOVTCR)  ----  cm  The  distance  in 

the  y-direction 
a  tracer 

moves  on  a  given 
subcycle  of 
INFACE. 


DMIN  =  Z  C 24)  -  -  INPUT  parameter. 

Allowable  relative 
error  in  energy 
sum.  When  rela¬ 
tive  error  is 
>  DMIN,  calcu¬ 
lation  is  termi¬ 
nated.  The  relative 
energy  error  is 
checked  in  EDIT. 
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Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


DNMS 

(UVCALC) 

-  -  -  - 

,  3 
g/cm 

The  average  den¬ 
sity  of  all  the 
"master"  materials 

in  a  slipline 
cell . 

DNSL 

(UVCALC) 

— 

/  3 
g/cm 

The  average  den- 

sity  of  all  the 
"slave"  materials 
in  a  slipline 
cell . 


DT  =  Z  (3) 


DTFACT  (SPHASE) 


DTMIN  =  Z (144)  ----  sec 


DTNA  =Z  (26)  -  sec 


DTSTR  (SPHASE)  -  sec 


DX  B.C.  (IMAX)  cm 


Time  step  At  = 
/ (AX.AY)min  \ 
\ (C+ (U,V)max) / 


•STAB. 

min 


Calculated  in  CDT. 


(See  CDT,  Section 

10.2) 


Factor  used  in 
calculating  a 
variable  time 
step  when  sub¬ 
cycling  the 
SPHASE  calculations. 


INPUT  parameter. 

Used  in  CDT  as  time 
step  cutoff.  After 
STAB  =  FINAL,  if 
DT  <  DTMIN  execution 
is  stopped. 

DT  from  previous 
time  cycle.  Used 
in  INPUT,  CDT,  EDIT, 
HELP  and  REZONE. 

The  time  step  used 
by  SPHASE:  DTSTR  = 
DT*DTFACT .  When 
SPHASE  is  not  sub- 
cycled,  DTFACT  =  1 
and  DTSTR  =  DT . 

The  radial- 
dimension  of  cells. 


Variable  Dimension 


Name 

Location 

of  Array 

Units 

Definition 

DXF 

=Z (136) 

cm 

DXF  =  DX(1)  if  the 
radial  dimension 
of  the  cells  is 
uniform. 

DXOO 

B.G. 

— 

cm 

DXOO  =  DX(0). 

DY 

B.C. 

( JMAX) 

cm 

The  axial- 
dimension  of  cells. 

DYF 

=Z (137) 

cm 

DYF  =  DY ( 1 )  if  the 
axial  dimension  of 
the  cells  is  uni¬ 
form. 

DYOO 

B.C. 

— 

cm 

DYOO  =  DY ( 0 )  . 

EAMMP 

t 

(T PHASE) 

ergs/g 

Specific  internal 
energy  of  mass 
transported  across 
the  right  edge  of 
a  cell. 

EAMPY 

(TPHASE) 

ergs/g 

Specific  internal 
energy  of  mass 
transported  across 
the  top  of  a  cell. 

ECK 

=Z(76) 

“  “  “  “ 

—  —  “  — 

Used  in  EDIT. 
Relative  error 
in  energy  sum: 

ECK=(£  Ek-ETH)/ETH 
k 

where  E^  is  total 
energy  in  cell  k. 

If  | ECK |  >  DMIN, 
execution  is 
stopped. 

EMIN 

=Z(85) 

ergs/g 

INPUT  parameter. 
Minimum  specific 
internal  energy 
to  be  used  in 
the  gamma -law 
equation  of  state. 
Usually  EMIN  *  107 

Variable  Dimension 

Name  Location  of  Array  Units  Definition 

EMIX  (ENCHCK)  -  ergs  The  sum  of  the 

energy  in  all  the 
cells  in  the  grid. 
Used  by  ENCHCK  to 
compute  the  rela¬ 
tive  energy  error 
after  each  phase 
of  the  cycle.  This 
energy  sum  is  com¬ 
puted  using  material 
velocities  rather 
than  cell-centered 
velocities  in 
interface  cells. 


EMOB 

-Z  (150) 

ergs 

Calculated  in  TPIIASE. 
Printed  in  EDIT. 

Total  energy  trans¬ 
ported  across  bot¬ 
tom  grid  boundary. 

EMOR 

®Z  (135) 

ergs 

Calculated  in  TPIIASE. 
Printed  in  EDIT. 

Total  energy  trans¬ 
ported  across  right 
grid  boundary. 

EMOT 

-Z  (146) 

ergs 

Calculated  in  TPIIASE. 
Printed  in  EDIT. 

Total  energy  trans¬ 
ported  across  top 
grid  boundary. 

ENERGY 

B.C. 

— 

ergs/g 

Defined  in  CDT  as 

the  specific  internal 
energy  of  a  cell  or 
of  a  material  in  a 
multimaterial  cell. 
Used  in  EQST  to 
compute  P  =  f (ENERGY, 
RIIOW)  . 

EOB  =  Z ( 1 3 4 )  -  ergs  Calculated  in  IIPHASE 

and  SPHASE.  Printed 
in  EDIT.  Change  in 
the  energy  of  the 
grid  due  to  work 
done  at  the  bottom 
boundary  of  the 
grid.  (See  Section 


Variable  Dimension 

Name  Location  of  Array  Units  Definition 

LOR  =Z(132)  -  ergs  Calculated  in  HPHASL 

and  SPHASE. 

Printed  in  EDIT. 
Change  in  the  energy 
of  the  grid  due  to 
work  done  at  the 
right  boundary  of 
the  grid.  (See 
Section  3.3.2) 

EOT  =  Z f 1 3 3 1  -  ergs  Calculated  in  HPHASE 

and  SPHASE. 

Printed  in  EDIT. 
Change  in  the  energy 
of  the  grid  due  to 
work  done  at  the  top 
boundary  of  the  grid 
(See  Section  3.3.2) 

ERDUMP  B.C.  -  -  Used  in  EDIT  and 

ERROR.  Flags  EDIT 
to  do  only  a  tape 
dump  on  an  error 
exit . 


ERR  (SPHASE)  (IMAX+1,3)  1/sec  The  normal  radial 

instantaneous 
strain  rate  associ¬ 
ated  with  cell 
centers.  (See 
Section  2. 3. 2.1) 


ERRINT  MKPLUG  ----  1/sec  The  normal  radial 

instantaneous 
strain  rate  cor¬ 
rected  for  con¬ 
vection.  Computed 
in  SPHASE,  used 
to  update  deviator 
stress,  Spj* ,  and  to 
compute  tne  angle 
of  maximum  shearing 
stress  in  PLGALF. 


Variable  Dimension 

Name  Location  of  Array 

ERZ  (SPHASE)  (IMAX+1,3) 


Units 


1/sec 


ERZINT 


MKPLUG 


1/sec 


(EQST) 


ergs/g 


(EQST) 


(30) 
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Definition 

The  instantaneous 
shear  strain  rate 
associated  with  cell 
centers.  (See 
Section  2 . 3 . 2 . 1) 

The  instantaneous 
shear  strain  rate 
corrected  for  con¬ 
vection.  Computed 
in  SPHASE,  used  to 
update  shear  stress, 
Srz,  and  to  compute 
the  angle  of  maxi¬ 
mum  shearing  stress 
in  PLGALF. 

Es,  the  specific 
internal  energy 
necessary  to  bring 
material  to  vapor 
temperature.  If 
material  is  ex¬ 
panded  and  its 
specific  internal 
energy  is  greater 
than  Es,  but  less 
than  Eg ,  an 
equation  of  state 
is  used  which  is  a 
blend  of  the  com¬ 
pressed  and  ex¬ 
panded  formulations. 
(See  Section  2. 3. 1.2) 

Defined  in  DATA 
Statement.  Values 
of  "a"  in  equation 
of  state  for  19  inert 
materials  fi  4  explosives, 
£  = ( y ~ 1 )  when  using 
perfect  gas  equa¬ 
tion  of  state. ] 


O  O,  4 


i  ^  ''J*  T  \ 


*9 


Variable  Dimension 

Name  Location  of  Array  Units  Definition 


ESALPH 

(EQST) 

(30) 

Defined  in  DATA 

Statement.  Values 
of  "a"  in  equation 
of  state  for  19 
inert  materials  and 

4  explosives. 

ESB 

(EQST) 

(30) 

Defined  in  DATA 

Statement.  Values 
of  "b"  in  equation 
of  state  for  19 
materials. 

ESBETA 

(EQST) 

(30) 

Defined  in  DATA 

Statement.  Values 
of  "3"  in  equation 
of  state  for  19 
inert  materials  and 

4  explosives. 

ESCAPA 

(EQST) 

(30) 

dynes/cm2 

Defined  in  DATA 

Statement.  Values 
of  MA"  in  equation 
of  state  for  19  inert 
materials  and  4 
explosives . 

ESC APB 

(EQST) 

(30) 

dyrtes/ cm2 

Defined  in  DATA 

Statement.  Values 
of  "B"  in  equation 
of  state  for  19  inert 
materials  and  4 
explosives . 

ESES 

(EQST) 

(30) 

ergs/g 

Defined  in  DATA 

Statement.  Values 
of  "E  "  in  equation 
of  stite  for  19  inert 
materials . 

ESESP 

(EQST) 

(30) 

ergs/g 

Defined  in  DATA 

Statement.  Values 
of  "E"'  in  equation 
of  st§te  for  19  inert 
materials . 
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Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


ESEZ 


ESP 


ESUM 


(EQST)  (30)  ergs/g  Defined  in  DATA 

Statement.  Values 
of  "E0"  in  equation 
of  state  for  19 
inert  materials. 

(EQST)  -  ergs/g  Eg,  the  specific 

internal  energy 
necessary  to 
vaporize  a 
material.  If  a 
material  is  ex¬ 
panded  and  its 
specific  internal 
energy  is  greater 
than  Eg  the  ex¬ 
panded  form  of 
the  Tillotson 
equation  of  state 
is  used.  (See 
Section  2. 3. 1.2) 

ergs  The  sum  of  the 

energy  in  all  the 
cells  in  the  grid. 
Used  by  EDIT  to 
compute  the  rela¬ 
tive  energy  error. 
In  ENCHCK,  ESUM 
is  the  energy 
computed  using 
cell  centered 
velocities  rather 
than  material 
velocities  in 
interface  cells. 
(See  EMIX.)  In 
UVMOD,  ESUM  is 
the  average 
specific  internal 
energy  of  a  slip¬ 
line  cell. 

p/p0,  the  com¬ 
pression  of  a 
material.  Used 
in  EQST  to  com¬ 
pute  the  material 
pressure . 


(EDIT, 

ENCHCK, 

UVMOD) 


ETA 


(EQST) 


vv  73.  r? 


»  st»  «* 


■^-g-^s'T  •^v^-';rT  "V^^Tr^T" 


•3777.V  *7777^77  .*.s,,.*.\»»T>l?,.’?l.l 
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Variable 

Name 

Location 

Dimens  ion 
of  Array 

Units 

Definition 

ETH 

=  Z  (13) 

ergs 

Theoretical  value 
of  total  energy  in 
the  mesh.  Calcu¬ 
lated  in  FILGRD 
initially. 

EVAPEN 

=  Z(101) 

— 

ergs 

Calculated  in 
TPHASE  and  RNDOFF. 

Printed  in  EDIT. 
Sum  of  energy  lost 
through  "evapo¬ 
ration"  of  mass 
left  in  cells  due 
to  round-off 
error . 


EVAPM  =  Z  (100)  -  g  Calculated  in 

TPHASE  and  RNDOFF . 

Printed  in  EDIT. 

Sum  of  mass  lost 
through  "evapo¬ 
ration".  See 
EVAPEN. 

EVAPMU  =  Z (102)  -  g-cm/sec  Calculated  in 

TPHASE  and  RNDOFF. 

Printed  in  EDIT. 

Sum  of  radial 
momenta  lost 
through  "evapo¬ 
ration"  .  See 
EVAPEN. 


EVAPMV  -Z (1 03)  -  g-cm/sec  Calculated  in 

TPHASE  and  RNDOFF. 
Sum  of  axial 
momenta  lost 
through  "evapo¬ 
ration".  See 
EVAPEN. 


EXPMIN 


(EQST) 


of 


EXPMIN- (p0)  is  the 
density  which  cor¬ 
responds  to  a 
minimum  value 
the  function 
Ap  +  Bp2  used 
the  Tillotson 
equation  of  state. 
In  order  for  the 


in 
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Variable 

Dimension 

Name 

Location 

of  Array 

Units 

Def init ion 

pressure  iteration 
to  converge,  P(p,E) 
must  have  a  positive 
slope.  (See 
Section  4 . 6 . 1 . 3) 

EZ  (EQST)  -  ergs/g  E0  in  the  thermal 

pressure  term  of 
the  Tillotson 
equation  of  state. 
(See  Section  2.  3. 1.2) 

EZPH2  =  Z  (104)  -  ergs  Sum  of  the  internal 

energy  increments 
that  are  set  to 
zero  in  TPHASE 
when  A  siE  < 

SIEMIN.  Printed 
in  EDIT. 

EZZ  (SPIIASE)  (IMAX+1,3)  1/sec  The  normal  axial  in¬ 

stantaneous  strain 
rate  associated  with 
cell  centers.  (See 
Section  2. 3. 2.1) 

EZZINT  MKPLUG  -  1/sec  The  normal  axial  in¬ 

stantaneous  strain 
rate  corrected  for 
convection.  Computed 
in  SPHASE,  used  to 
update  deviator 
stress,  S  ,  and  to 
compute  the  angle 
of  maximum  shearing 
stress  in  PLGALF. 

FB  (DMADJ)  -  g  The  mass  being 

transported  out 
of  a  cell  across 
its  bottom  bound¬ 
ary.  Computed  and 
used  by  DMADJ  to 
exactly  evacuate 
materials  from 
interface  cells. 
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Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

FINAL 

=  Z (113) 

INPUT  parameter. 
Maximum  value  of 
stability  fraction 
(STAB).  If 

FINAL  =  0,  the 
stability  fraction 
will  be  constant. 
Used  in  CDT. 

FL 

(DMADJ) 

g 

The  mass  being 
transported  out 
of  a  cell  across 
its  left  boundary. 
(See  FB) 

FLEFT 

6 

B.C. 

(JMAX) 

g-cm/sec 

Radial  momentum  of 
mass  transported 
across  left  side 
of  cell.  Equiva- 
lenced  to  UL  array. 
(See  Section 

10.4. 2) 

FLX 

(DMADJ) 

(4) 

g 

The  mass  trans¬ 
port  terms  for  a 
given  material  at 
the  four  boundaries 
of  a  cell.  Com¬ 
puted  and  used  by 
DMADJ  to  adjust 
transport  terms 
so  that  materials 
are  not  over¬ 
emptied  . 

FR 

(DMADJ) 

g 

The  mass  being 
transported  out 
of  a  cell  across 
its  right  bound¬ 
ary.  (See  FB) 

FRACRT 

B.C. 

(NVOID, 

NMXCLS) 
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n 

cm"' 

The  fractional 
area  of  the  right 
cell  boundary 
associated  with  a 
given  material  in 
an  interface  cell. 
Used  to  compute 
mass  transport 
terms  of  that 
material . 

EVfVift  “% 


Variable 

Name 

FRACTP 


Location 


B.C. 


Dimension 
of  Array 

(NVOID, 

NMXCLS) 


FRACX 


B.C. 


FRACY 


(FRACS) 


(MOVTCR, 
PLGGEN , 
PLUGUV, 
FRACS) 


(MOVTCR, 

PLGGEN, 

PLUGUV, 

FRACS) 


Units  Definition 

2 

m  The  fractional 

area  of  the  top 
cell  boundary 
associated  with  a 
given  material  in 
an  interface  cell. 
Used  to  compute 
mass  transport 
terms  of  that 
material . 

—  The  fraction  of 
the  x-dimension  of 
the  cell  from  the 
left  of  the  cell 

to  the  intersection 
of  the  interface 
with  the  cell's  top 
boundary.  Used  in 
NEWRHO,  as  well  as 
in  FRACS,  to  in¬ 
dicate  which 
boundary  of  a  cell 
is  crossed  by  the 
interface.  (Equiv- 
alenced  to  1VSC  in 
FRACS.) 

The  fraction  of 
the  y-dimension 
of  the  cell  from 
the  bottom  of  the 
cell  to  the  inter¬ 
section  of  the  inter 
face  with  the 
cell's  right 
boundary. 

---  The  fractional 

part  of  the  x- 
coordinate  of  a 
tracer  (in  cell 
units). 

—  The  fractional 
part  of  the  y- 
coordinate  of  a 
tracer  (in  cell 
units) 
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Variable 

Dimension 

■  Name 

Location 

of  Array 

Units 

Definition 

s 

(DMADJ) 

— 

g 

The  mass  being 
transported  out 
of  a  cell  across 

GAMC 


B.  C. 


(JMAX) 


GAMGAS 


(CDT) 


GAMMA 


;Z  (62) 


HOOP 


(SPHASE) 


dynes/cm' 


=Z (88) 


its  top  boundary. 
(See  FB) 

Mass  transported 
across  the  left  side 
of  a  cell.  Equiva- 
lenced  to  PL 
array.  (See 
Section  10.4.2) 

A  constant  used 
in  approximating 
the  sound  speed 
of  a  high  explo¬ 
sive.  Defined  in 
a  DATA  statement 
in  CDT  for  explo¬ 
sives.  For  ideal 
gas  CDT  sets 
GAMGAS ( 20) =GAMMA. 

INPUT  parameter. 

"y"  in  the  ideal 
gas  equation  of 
state: 

P  =  (y-I)pE.  Used 
in  EQST  and  CDT. 

s  ee  T  ’  (Srr+Szz)  * 
the  hoop  stress, 

calculated  when 
in  cylindrical 
coordinates . 

Used  in  most 
subroutines  as  the 
index  on  the  grid 
columns . 

A  calling  argument 
in  DMCALC  which 
identifies  the 
leftmost  column 
of  the  grid  for 
which  the  mass 
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teMKISSS 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


ICLADD  =Z (38) 


ICSTOP  =Z  (7)  - 

ICY  INF  - 

IDL  (MAP, TPHASEJ  ---- 


transport  terms  of 
interface  cells 
are  to  be  computed 
(DMCALC  is  called 
from  NEWMIX  for  a 
single  cell  and 
from  INFACE  for 
all  interface  cells 
in  the  active  grid.) 

Flag  used  when 
detonating  a  high 
explosive.  When 
ICLADD  =  -  1 
either  all  explo¬ 
sives  have  been 
detonated  or  the 
problem  does  not 
involve  explosives. 
(Defined  in  DETIME 
initially,  later 
in  ADDENG.)  When 
ICLADD  =  0  some 
explosive  in  the 
problem  has  not 
yet  been  detonated. 

INPUT  parameter. 

Used  in  EDIT. 

Execution  stops  on  the 
ICSTOP  cycle  when 
stopping  on  a  speci¬ 
fied  cycle  rather 
than  on  time. 

The  number  of  passes 
through  INFACE 
(= INT (CYCMX) ) . 

The  rightmost  column 
of  the  grid  repre¬ 
sented  by  the  sym¬ 
bolic  maps .  On 
cycle  0,  I DL=  JMAX ; 
thereafter  IDL=I1, 
the  radial  active 
grid  counter. 


-  *-  »  «-  »  *— «  Ifc-  fc  -m  *■■*.  »«*  ■»*  - 


fcT  ifci  i^a.i 


Variable 

Name 


Location 


Dimension 
of  Array 


Units 


I DLT 


m  * 


-Z(19) 


I DRT 


=  Z  (  3 1 ) 


IEXTX 


=  Z  (123) 


I  FLUX 


(DMCALC) 
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Definition 


The  leftmost  column 
of  the  grid  that 
contains  a  high 
explosive.  Defined 
in  DETIME. 


The  rightmost  column 
of  the  grid  that 
contains  a  high 
explosive.  Defined 
in  DETIME.  When 
I DRT  =  0,  the 
calculation  does 
not  involve  an  ex¬ 
plosive  . 


A  calling  argument 
in  DMCALC  which 
indentifies  the 
right-most  column 
of  the  grid  for 
DMCALC  to  compute 
(see  IB) . 


Used  in  REZONE  and 
TPHASE.  When 
IEXTX  =  1  the  grid 
is  rezoned  in  the 
x-direction.  Enables 
user  to  rezone  in 
the  x-direction  only. 
Automatic  rezones 
will  not  be  trig¬ 
gered  by  signals 
reaching  the  right 
grid  boundary  unless 
IEXTX  =  1. 


In  DMCALC,  when  IFLUX= 
1,  the  top  and  right 
mass  transport  terms 
are  computed;  when 
IFLUX=2,  the  bottom 
and  left  terms  are 
computed  for  the  cell 
in  column  IT  and  row 
JT. 


Variable  Dimension 

Name  Location  of  Array  Units  Definition 

IFS1  (DMCALC)  -  -  When  IFS1=1,  the 

cell  (IT, JT)  DMCALC 
is  operating  on  is 
a  free  surface  cell. 
This  cell  will  not 
be  the  "donor"  cell 
unless  the  neighbor 
cell  is  also  a  free 
surface  cell  (IFS2=1). 

I FS2  (DMCALC)  -  -  When  IFS2  =  1,  the 

neighbor  cell  DMCALC 
is  operating  on  is  a 
free  surface  cell.  This 
neighbor  cell  will  not 
be  the  "donor"  unless 
the  cell  (IT,JT)  is 
also  a  free  surface 
cell  (IFS1-1). 

IGM  =  Z  (21)  -  -  A  flag  which  indicates 

which  coordinate  sys¬ 
tem  is  being  used: 

IGM=1  for  plane 

coordinates 
IGM=0  for  cylindrical 
coordinates. 

IMAX  =Z(33)  -  -  INPUT  parameter. 

Number  of  columns 
in  the  mesh.  IMAX 
must  be  an  even 
number  if  the  grid 
is  to  be  rezoned  in 
the  x-direction. 

INTER  =  Z  ( 8 7 )  -  -  A  flag  which  causes 

certain  diagnostic 
messages  to  be 
printed  each  cycle. 

(See  Section  9.2) 

IPCYCL  =Z(49)  -  -  INPUT  parameter.  Used 

in  EDIT.  The  number 
of  cycles  between 
EDIT  prints  when  edit¬ 
ing  on  cycles  rather 
than  on  time. 
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Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

IPLGBT 

•Z(60) 

The  bottom-most  row 
of  the  plugging 
region  of  the  target. 
Used  only  in  plugging 
calculations . 

IPLGRT 

=Z (56) 

The  right-most  row 
of  the  plugging 
region  of  the  target. 
Used  only  in  plugging 
calculations . 

IPLGTP 

=  Z (61) 

The  top-most  row 
of  the  plugging 
region  of  the  target. 
Used  only  in  plugging 
calculations . 

IPR 

= z  CIS) 

Maximum  number  of 
iterations  to  be  per¬ 
formed  by  CDT  to 
achieve  pressure 
equilibration  between 
materials  in  multi¬ 
material  cells. 

ITC 

(CDT) 

— 

— 

The  number  of  itera- 

X  ✓  . 

tions  CDT  has  processed 
while  equilibrating 
the  pressure  of 
materials  in  a  multi¬ 
material  cell.  If 
ITC  exceeds  the  input 
cutoff,  IPR,  the  cal¬ 
culation  is  stopped. 

ITFLAG  (ADDTCR)  -  -  A  flag  used  in  ADDTCR. 

ITFLAG=0  when  points 
are  interpolated  in 
physical  coordinates 
(NADD>0) 

ITFLAG=1  when  points 
are  interpolated  in 
cell  units  (NADD<0), 

The  second  option  is 

used  only  when  the  cell 
dimensions  are  constant, 
in  which  case  SLTUPA  sets 
NADD  negative. 
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Variable 

Name 

ITXB 


ITXl 


Location 


(FRAGS) 


(FRAGS) 


I)  i  mens  ion 
of  Array 


Units 


ITX2 


(FRAGS) 


ITYB 


ITY1 


(FRAGS) 


(FRAGS) 


ITY2 


IVARDX 


IVARDY 


(FRAGS) 


-Z  (83) 


=Z(54) 


=  Z  (47) 
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Definition 

ITXB+1 is  the  column 
that  contains  the  first 
tracer  of  a  pair  being 
considered  by  FRAGS. 

The  vertical  grid  line 
crossed  by  the  inter¬ 
face  for  which  FRACS 
is  computing  a  frac¬ 
tional  cell  face  area. 

ITX2+1  is  the  column 
that  contains  the 
second  tracer  of  a 
pair  being  considered 
by  FRACS. 

ITYB+1  is  the  row  that 
contains  the  first 
tracer  of  a  pair  being 
considered  by  FRAGS. 

The  horizontal  grid 
line  crossed  by  the 
interface  for  which 
FRACS  is  computing  a 
fractional  cell  face 
area . 

ITY2+1  is  the  row  that 
contains  the  second 
tracer  of  a  pair  being 
considered  by  FRACS. 

When  1 VARDX= 1 ,  the 
radial  dimension  of 
cells  is  non-uniform. 

When  IVARDY=1,  the 
axial  dimension  of 
cells  is  non-uniform. 

INPUT  parameter.  II 
is  used  to  limit  the 
calculation  in  the 
radial  direction  to 
an  "active  mesh". 
Initially,  11=2  + 


is 


■<i 


Variable  Dimension 

Name  Location  of  Array 


Units 


Definition 


12 


13 


=Z(48) 


B.C. 


=Z (89) 


JB 


JDBT 


=Z (28) 


JDL  (MAP,TPHASE) 
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the  column-number 
of  the  last  column 
in  which  a  velocity 
and/or  internal 
energy  exists.  Up¬ 
dated  in  HPHASE, 
TPHASE,  ADDENG,  and 
REZONE. 


INPUT  parameter. 
Like  II  but  for  the 
axial  direction. 


The  rightmost  column 
of  cells  printed  by 
EDIT.  13=1  for  a 
"short"  print,  13=11 
(the  radial  edge  of 
the  active  grid)  for 
a  "long"  print. 

Used  in  most  sub¬ 
routines  as  index 
on  the  grid  rows. 

A  calling  argument  in 
DMCALC  which  identifies 
the  bottom-most  row 
of  the  grid  for  which 
the  mass  transport 
terms  of  interface 
cells  are  computed. 
(DMCALC  is  called  from 
NEWMIX  for  a  single 
cell  and  from  INFACE 
for  all  interface  cells 
in  the  active  grid.) 

The  bottom  row  of  the 
grid  that  contains  a 
high  explosive.  De¬ 
fined  in  DETIME. 

The  topmost  row  of  the 
grid  represented  by 
the  symbolic  maps.  On 
cycle  0,  JDL=JMAX; 
thereafter,  JDL=I2, 
the  axial  active  grid 
counter . 


Variable 

Dimension 

Name 

Location 

of  Array 

Units 

De  f ini t ion 

JDTP 


=  Z  (29) 


The  topmost  row  of  the 
grid  that  contains  a 
high  explosive.  De¬ 
fined  in  UETIME. 

A  calling  argument  in 
DMCALC  which  identifies 
the  topmost  row  of  the 
grid  for  DMCALC  to 
compute  (see  JB) . 


JEXTY 


=  Z  (124) 


A  flag  used  in 
REZONE.  When 
JEXTY  =  1,  the 
grid  is  rezoned 
in  the  y- 

direction.  Enables 
user  to  rezone  in 
y-direction  only. 
Automatic  rezones 
will  not  be  trig¬ 
gered  by  signals 
reaching  the  upper 
grid  boundary 
unless  JEXTY  =  1. 


JMAX 


=  Z  (35) 


JSLIP 


(PLGTCR) 


INPUT  parameter. 
Number  of  rows  in 
the  grid.  JMAX 
must  be  an  even 
number  if  the  grid 
is  to  be  rezoned 
in  the  y-direction. 

The  row  across  which 
the  plug  is  being 
extended. 


=  Z (90) 


B.C. 


Used  as  cell-index  in 
all  subroutines. 

K=  (J- 1)  MMAX+1  +  1. 

Used  as  cell-index 
for  the  cell  above 
K.  KA  =  K  +  IMAX. 
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Variable  Dimension 

Name  Location  of  Array 

KB  (DM CALC , SPHASE)  ---- 


(EDIT, INPUT, 
SETUP) 


KMAX 


=  Z  (37) 


(EDIT, INPUT, 
SETUP) 


KSPACE 


(EDIT) 


(EDIT, INPUT, 
SETUP) 


KUNITR 


=Z(14) 


Units  Definition 

-  Used  as  cell-index 

for  the  cell  below 
cell  K.  KB-K-IMAX. 

-  The  number  of  words 

in  the  DETIM  array 
written  on  the  re¬ 
start  file*  KDT  =  1 
if  there  is  no  ex¬ 
plosive  in  the 
problem;  otherwise, 

KDT  =  KMAX. 

-  The  K- indices  of 

the  four  overlap 
cells  used  to 
compute  the  velocity 
components  of  a 
tracer  particle. 

----  Calculated  in  SETUP 

(=IMAX*JMAX+1) .  K- 
index  of  the  last 
cell  in  the  grid. 

-  The  number  of  words  in  the 

PLWP  array  written  on  the 
restart  file.  KPW=1,  if 
the  plugging  option  is  not 
being  used;  otherwise, 
KPW=NTCC,  the  number  of 
passive  tracers. 

-  Used  for  spacing  printed 

output . 

-  The  number  of  words  in 

the  STRSZZ,  STRSRR  and 
the  STRSRZ  arrays  written 
on  the  restart  file. 

KSS=1  if  the  strength 
phase  is  omitted  (CYCPH3* 
-1);  otherwise ,KSS=KMAX. 

_  Number  of  the  file 

read  by  INPUT. 
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V 

V..'* 


M 

m 


z ’r*y«v'*yrir*.  -^"-’7 


"7?7>: 


p,v  »v 


.■rij.lyyii  Ty»y  ■»rrnf  yypryy  M-  *J»  *.">i " 


Variable 

Name 


In  mens  ion 

Location  of  Array 


Uni  ts 


Definition 


KUNITW 


=Z (17) 


LABLH 


=  Z(1) 


U) 


LAST 


B.C, 


LJ 


INF 


LSAVE 


(CUT , PLGTCR) 


LVF 


(MOVTCR, 

ENDMV) 


LVISC 


:  Z  (116) 


M 


=Z (9 1) 
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Number  of  the  file 
written  on  by  SETUP 
and  EDIT . 


An  array  used  by  CARDS 
to  store  integer  variables 
in  the  Z-block  of  blank 
common . 


The  last  word  of  blank 
common.  Used  in  the  main 
routine  to  initialize 
the  blank  common  storage 
to  zero. 


The  number  of  the  INFACE 
subcycle.  (LJ  goes  from 
1  to  ICY.) 


In  CDT  -  the  material 
package  number  of  the  only 
material  in  an  interface 
cell  that  contains  only 
one  material. 


In  PLGTCR  -  the  second 
index  of  the  target  tracer 
to  the  right  of  the  inter¬ 
cept  of  the  plug  edge  with 
the  back  surface  of  the 
target . 


The  second  index 
void  tracer  that 
ponds  to  a  given 
tracer . 


of  a 
cor res  - 
material 


A  linear  artificial  vis¬ 
cosity  term  is  added  to 
the  cell  boundary  pres¬ 
sures  when  LVISC=1. 


Usually  the  index  which 
links  the  material  arrays 
to  cell  K: 


M  =  MFLAG ( k) - 100  and 
AMX(K)=£XMASS(N,M). 


Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

MA 

B.C. 

— 

— 

Usually  the  value 
of  MFLAG  for  the 
cell  above  cell  K. 

MAPS 

=  Z(42) 

— 

— 

INPUT  parameter 
that  determines  the 

printing  of  symbolic 
maps  on  EDIT  cycles. 

If  MAPS  =  0,  maps 
are  not  printed;  if 
MAPS  =  1,  they  are 
printed  and  the 
first  map  is  of 
compression;  if 
MAPS  =  2,  they  are 
printed  and  the 
first  map  is  of 
density . 

MASTER  (FRACS)  -  -  While  processing 

the  package  n  inter¬ 
face,  FRACS  sets 
MASTER  B  n  if  package 
n  is  a  "master" 
package.  Its  tracers 
will  be  used  to 
compute  the  angle  of 
the  slipline  in  each 
slipline  cell. 

MASTRD  SLIDPK  (NMAT)  -  If  MASTRD(n)  =  n, 

then  package  n  is  a 
master  package  and  all 
or  part  of  the  package 
n  interface  is  a  slip¬ 
line.  If  MASTRD (n) =0 , 
then  package  n  is  a 
slave  package 
(NSLAVD(n)  =  n)  or 
it  is  neither  ancPno 
part  of  its  interface 
is  a  slipline. 


V 


r  • 

3 


Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


MAT 


MXCELL 


(30) 


MAXX 

MAXY 


=  Z  (120) 
*Z  (122) 


MFLAG 


B.C. 


( KMAX) 


MINX 

MINY 


=Z (119) 
=  Z  (121) 


MO 


B.C. 
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Material  code  numbers 
for  the  packages. 
MAT(1)=3  indicates 
the  material  in 
package  one  is  iron. 
See  list  in  Section 
2.3.1  of  this  report 
or  comments  in  the 
EQST  subroutine. 


f  Used  in  ADDTCR.  Right 
column  and  top  row, 
respectively,  of  the 
region  in  which  tracer 
oarticles  are  to  be 
added  (see  NADD,  MINX, 
MINY) . 


A  flag  which  indicates 
whether  a  cell  is  pure 
or  contains  an  inter¬ 
face.  If  MFLAG (K)  < 
100,  cell  K  is  pure. 

If  MFLAG (K)  >  100, 
cell  K  contains  an 
interface  and  the 
flag  gives  the  storage 
location  in  the 
material  arrays  for 
cell  K. 


fused  in  ADDTCR.  Left 
column  and  bottom 
row,  respectively,  of 
the  region  in  which 
tracer  particles  are 
to  be  added  (see  NADD, 
l MAXX,  MAXY). 


The  MFLAG (K)  of  cell 
K  when  it  was  pure. 
Used  by  NEWMIX  to 
define  the  material 
variables  (XMASS,  SIE, 
US,  VS,  etc.)  of  a 
cell  that  has  just  be¬ 
come  an  interface  cell 


awyMBia  i*  i  t  t  f  ^  t*-*-*  -.  ^  ^  l%  ■>  ^ 


Var iab le 
Name 


Dimension 

Location  of  Array 


Units 


Definition 


MOS 


B.C, 


MPLUS 


(Al)DTCR) 


MR 


B.C. 


MSLD 


SLIDE 


(NSLD) 


MTK 


MOV 


(4) 


N 


=  Z  (92) 


NAD 


(ADDTCR, 

EDIT.PLGADD) 
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The  number  of  cells  con¬ 
taining  the  slipline  on 
any  one  cycle.  Incre 
mented  in  PTSAV  which 
saves  the  slipline  inter¬ 
cepts  with  cell  boundaries. 


2 


An  accumulative  number 
of  tracers  added  to  a 
subpackage  or  package. 


Usuall>  the  value  of 
MFLAG  for  the  cell  on 
the  right  of  cell  K. 


A  flag  word  for  each 
slipline  cell  which 
stores  its  I,J,  and 
MFLAG (K)  values: 


MSLD (MS) “MFLAG (K) *10* 

+  ino3  +  j  . 


This  allows  the  user 
to  have  smaller  di¬ 
mensions  for  the  slip¬ 
line  arrays XINT,  YINT 
than  for  the  THETA 
array  which  is 
associated  with 
every  interface  cell. 


The  MFLAG  values  of 
the  four  overlap  cells 
used  to  compute  the 
velocity  components 
of  a  tracer  particle. 


In  EQST,  N  is  the 
material  code  number 
transferred  from  CDT. 


Temporary  storage  for 
NADD  in  ADDTCR.  In 
PLGADD,  NAD  is  the 
number  of  tracers 
added  at  the  corner 
where  the  plug  will 
begin. 

I NADD I 


In  EDIT,  NAD 


1 


•1 


a 


* 


Variable 

Name 

Location 

Dimens i on 
of  Array 

Units  Definition 

NADD 

=  Z  ( 1 1 8) 

-  NADD  has  two  functions. 

In  EDIT  it  controls  the 
frequency  of  calls  to 
ADDTCR.  In  ADDTCR  it  is 
used  to  determine  how 
many  tracers  to  add  be¬ 
tween  two  existing 
tracers.  Its  input  value 
is  saved  in  NAD  and  re¬ 
stored  at  the  end  of 
ADDTCR. (See  Section  8.2) 

NAME 

(ENCHCK) 

-  A  string  of  hollerith 

symbols  passed  by  HELP 
through  a  calling  argu¬ 
ment  which  identifies 
which  phase  of  the  calcu¬ 
lation  (SPHASE ,HPHASE , or 
TPHASE)  has  just  been 
completed. 

NBGMD 

SLIDPK 

(NMAT) 

The  second  index  of  the 
tracer  which  begins  the 
slipline  of  a  master 
package . 

NBGP 

(THETAS) 

----  The  slave  package  whose 

"beginning"  slipline 
tracer  (NBGSD)  is  an 
endpoint  of  the  slipline. 

NBGSD 

SLIDPK 

(NMAT) 

-  The  second  index  of  the 

tracer  which  begins  the 
slipline  of  a  slave 
package . 

NC 

-1 (30) 

----  Cycle  number.  Set  ini¬ 

tially  to  -1  in  INPUT. 
Incremented  thereafter 
in  CDT.  (Decremented  in 

INPUT  when  restarting  a 
calculation . ) 

NDUMP7 

=  Z(6) 

— 

-  INPUT  parameter.  Used 

in  LDIT  to  control  fre- 

quency  of  restart  dumps, 
e.g.,  a  dump  will  occur 
on  every  LDIT  pi int  when 
NDUMP7  =  1 ,  or  on  eveTy 
fifth  EDIT  print  when 
NDUMP7  =  5. 
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Variable 

Name 

Location 

imens ion 
;.»f  Array 

Units  Definition 

NECYCL 

=  Z (77) 

-  Defined  and  printed  in 

EDIT.  The  cycle  on  which 
the  largest  relative 
error  in  the  energy  sum 
was  computed. 

NEND 

(ENDMV) 

_  -  -  - 

-  The  second  index  of  the 

plug  tracer  which  is  at 
the  top  of  the  slipline. 

NENDMD 

SLIDPK 

(NMAT) 

-  The  second  index  of  the 

tracer  which  ends  the 
slipline  of  a  master 
package. 

NENDMI 

l 

=  Z  C 1 0  8 ) 

-  The  second  index  of  what 

is  initially  the  last 
target  tracer  defining 
the  si ip line.  The  slip- 
line  gets  extended  along 
the  target  boundary  to 
allow  the  projectile  lip 
to  slip  on  the  target. 
Used  only  in  plugging 
calculation ; . 

NENDSD 

SLIDPK 

(NMAT) 

-  The  second  index  of  the 

tracer  which  ends  the 
slipline  of  a  slave 
package . 

NENP 

(THETAS) 

-  The  slave  package  whose 

"ending"  slipline  tracer 
(NENDSD)  is  an  endpoint 
of  the  slipline. 

NERR 

B.C. 

-  Used  in  ERROR  as  exit 

flag.  Prevents  ERROR 
from  being  called  more 
than  once  during  a 
single  run. 

NFLAG 

(CDT) 

— 

-  A  counter  used  in  the 

pressure  iteration  to 
count  the  number  of 
times  the  weighting 
factor,  CSQR,  is  doubled 
to  give  all  materials 


Variab le 
Name 

Location 

Dimension 
of  Array 

Units 

Definition 

a  positive  partial 
volume.  If  NFLAG 
exceeds  the  input  cut¬ 
off,  IPR,  the  calcula¬ 
tion  is  stopped. 

NFRELP 

=  Z(5) 

INPUT  parameter.  Used 
in  EDIT  to  control  fre¬ 
quency  of  "long"  prints. 

A  "long"  print  will 
occur  with,  every 

EDIT  if  NFRELP  =  1; 
when  every  fifth  EDIT 
if  NFRELP  =  5. 

NK 

B.C. 

“  “  “  “ 

-  -  -  - 

Tells  which  statement 
number  of  a  subroutine  is 
near  the  call  to  ERROR. 

NLINER 

=Z (1 05) 

The  package  number  of  the 
material  that  forms  the 
jet  of  a  shaped  charge. 
Used  only  when  calculat¬ 
ing  the  collapse  of  a 
shaped  charge  liner. 

NMAT 

-Z (68) 

The  number  of  material 
packages  in  the  problem. 
(Note:  each  material 

package  can  be  made  up 
of  several  disconnected 
subpackages.  Also 
two  packages  can 
contain  the  same  material 
but  have  different  ini¬ 
tial  conditions.) 

NMP 

B.C. 

(NMAT) 

The  number  of  tracer 
particles  used  to  circum¬ 
scribe  the  boundary  of  a 
given  material  package. 

NMXCLS 

-Z  (73)  . 

— 

— 

The  maximum  number  of 

interface  cells  that 
SETUP. or  NHWMIX  can 


generate.  This  number 
should  coincide  with  the 
dimension  of  variables  in 
the  MXCELL  common  block. 


rxr.?:ry7'f"^s.ry?$  ■’>  '  >  y^vx  y 


■  ’  TT-.'r^,^v  1 «!^7T VH7: yr^7 T v'A^; 


Variab le 
Name 


Locat i on 


Dimension 
of  Array 


Units 


Def ini tion 


NODUMP 


=Z(96) 


NOSLIP 


=  Z (115) 


NPRINT 


B.C. 


NR 


B.C. 


NRC 


(TPHASE) 


NREC 


(INPUT) 


NRT 


(TPHASE) 


NSAVE 


(FLGSET) 


NSLAVD 


SLIDPK 


(NMAT) 
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INPUT  parameter.  Used 
in  EDIT.  When  NODUMP=l 
no  restart  dumps  are 
written  except  by  SETUP 
on  cycle  0.  Allows  user  to 
restart  a  problem  with¬ 
out  writing  on  the 
restart  file. 


A  flag  which  when  set 
to  1  indicates  the 
calculation  has  no 
s lip  lines . 


NPRINT=1  during  the 
cycle  on  which  EDIT 
prints  and  checks  the 
energy  error. 


Identifies  which  sub¬ 
routine  calls  ERROR. 
Used  in  ERROR  for 
printing  error  message. 


Used  in  TPHASE  to  signal 
that  the  active  grid  must 
be  enlarged  one  column 
on  the  right. 


The  number  of  logical 
records  in  each  re¬ 
start  dump. 


Used  in  TPHASE  to  signal 
that  the  active  grid 
must  be  enlarged  one 
row  at  the  top. 


The  package  number  of  a 
material  in  a  cell  which 
was  generated  as  an 
interface  cell  but  is 
flagged  a  pure  cell  by 
FRACS  and  FLGSET. 


If  NSLAVD (n) =n ,  then 
package  n  is  a  slave 
package  and  all  or  part 


Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

of  the  package  n 
interface  is  a  slip¬ 
line.  If  NSLAVD (n) * 

0,  then  package  n 
is  either  a  master 
package  (MASTRD (n) = 
n)  or  no  part  of 
its  interface  is  a 
slipline. 

NSLD 

=  Z  (112) 

The  maximum  number 
of  slipline  cells 
in  the  grid  on  any 
one  cycle.  This 
number  should  coin¬ 
cide  with  the  di¬ 
mension  of  the 
variables  in  the 

SLIDE  common  block. 

NSLPMN 

=Z(107) 

The  initial  value  of 

NBGSD (2) .  (Used  only 
in  plugging  calcula¬ 
tions)  .  Updated  in 

PLGTCR  so  that  NSLPMX- 
NSLPMN+1  is  the  number 
of  identical  tracers  at 
the  corner  of  the  plug, 
as  the  plug  evolves. 

NSLPMX 

=  Z  (106) 

The  initial  value  of 
NENDSD(2) .  Used  in 
plugging  calculations 
only . 

NTCC 

=  Z  (81) 

INPUT  parameter.  When 
NTCC>0  passive  cell 
centered  tracers  are 
used. 

NTPMX 

=Z (78) 

— 

— 

INPUT  parameter.  Equals  . 

the  maximum  number  of 
tracers  that  SETUP  or 
ADDTCR  can  generate  to 
circumscribe  a  single 
material  package.  The 
value  of  this  number 
should  coincide  with  the 
dimension  of  the  TX  and 
TY  arrays. 
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Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


NTRACR  =  Z  (72)  -  -  INPUT  parameter. 

Used  to  determine 
the  density  of 
material  tracer 
particles.  (See 
Section  8.2) 

NUMREZ  =Z (12)  -  -  INPUT  parameter. 

Initially  equals 
the  number  of  auto¬ 
matic  rezones 
allowed  during  a 
calculation.  Dimin¬ 
ished  by  one  after 
each  rezone. 

NUMSCA  =Z(43)  -  -  INPUT  parameter.  Number 

of  times  the  print  inter¬ 
val  is  to  be  rescaled. 
Used  in  EDIT.  See  PRDELT 
for  further  details. 

NUMSP  -Z(4)  -  -  Used  in  EDIT  to  count 

the  number  of  "short" 
prints  since  the  last 
"long"  print.  When 
NUMSP-NFRELP,  EDIT  pro¬ 
duces  a  "long"  print, 
i*e.,  information  on 

all  the  cells  in  the 
active  grid,  rather  than 

just  on  the  cells  in  the 
first  column. 

NUMSP  T  =  Z  (41)  -  -  Used  in  EDIT  to  count 

the  number  of  prints 
(short  or  long)  since 
the  last  tape  dump. 

NUNxxx  =Z(xxx)  -  -  Unused  Z-storage. 

NVOID  B.C.  -  -  Defined  in  INPUT:  =NMAT 

+  1,  number  of  material 
packages  +1.  The  number 
of  the  void  package  if 
there  is  one. 
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Variab le 
Name 


Location 


Dimension 
of  Array 


Units 


Definition 


NVP 


(MOVTCR) 


NVRTEX 


=  Z (109) 


N10 


(CDT) 


Nil 


(CDT) 


OKE 


(IIPHASE) 


ergs 


OUTFLX 


(DMCALC) 


g 
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The  number  of  material 
packages  plus  one.  Used 
as  a  limit  on  the  do- loop 
which  moves  the  tracer 
particles.  The  void 
tracers  are  moved  simul¬ 
taneously  with  the 
material  tracers  and  the 
passive  tracers  are 
moved  (on  the  last  sub¬ 
cycle  of  INFACE)  when 
N=NVP. 


The  second  index  of  the 
void  tracer  at  the  vertex 
of  the  void  closing  re¬ 
gion.  If  NVRTEX=0 ,  the 
code  does  not  call 
VDCLOS  to  automatically 
close  a  void. 


Used  in  CDT  to  identify 
the  I-index  of  the  cell 
which  controls  the 
time  step. 


Used  in  CDT  to  identify 
the  J-index  of  the  cell 
which  controls  the 
time  step. 


The  kinetic  energy  of 
a  slipline  cell  before 
its  material  velocities 
are  updated  in  UVCALC. 


The  total  mass  of  a 
material  being  trans¬ 
ported  out  of  an  in¬ 
terface  cell.  Computed 
by  DMADJ  to  see  if  the 
mass  going  out  exceeds 
the  mass  in  the  cell. 
The  mass  transport 
terms  are  adjusted 
accordingly  to  prevent 
negative  masses. 


Variable 

N  ame 

Location 

Dimension 
of  Array 

Units 

! 

Definition 

OUTIE 

B.C. 

(NMAT) 

ergs 

The  total  internal  energy 
lost  from  each  package 
due  to  "evaporation", 
cut-offs,  or  material 
leaving  the  grid. 

Printed  by  EDIT. 

OUTKE 

B.C. 

(NMAT) 

ergs 

The  total  kinetic  energy 
lost  from  each  package 
due  to  "evaporation", 
cut-offs,  or  material 
leaving  the  grid. 

Printed  by  EDIT. 

P 

80 

n 

(KMAX) 

dynes /cm2 

1 

Cell  pressure.  Note 

that  the  P-storage 
is  used  by  other  arrays 
in  INFACE,  REZONE,  and 
TPHASE .  See  equi¬ 
valence  statements  in 
the  "included"  element. 

PABOVE  (HPHASE)  ----  dynes/cm2  The  pressure  at  the 

boundary  between  cell  k 
and  the  cell  above  it, 
an  average  of  the  cell- 
centered  pressures  of 
the  two  cells. 


PAY 


(CDT) 


PBLO 


(HPHASE) 


PI 

PIDY 


-Z(8) 
=  Z(8) 


dynes/cm' 


The  weighted  average 
pressure  computed  in  the 
pressure  iteration  to 
which  the  pressures  of 
all  materials  must 
converge.  (See  Section 
4. 6. 1.2) 


dynes/cm" 


The  pressure  at  the 
boundary  between  cell  k 
and  the  cell  below  it, 
an  average  of  the  cell- 
centered  pressures  of 
the  two  cells. 
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am 


z 


« 


Var iab le 
Name 


Location 


Dimension 
of  Array 


Units 


B.C. 


B.C. 


(See  dynes/cm' 

definition 
of  JDX2 
parameter 
in 

Appendix 

A) 


PLGOPT 


=  Z  (111) 


PLUGON 


=  Z  (125) 


MXCELL 


(NMAT) 


ergs 
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Definition 


The  first  array  follow¬ 
ing  the  150  word  Z-block 
PK(1)  =  problem  number, 
PK(2)  =  restart  cycle, 
PK(3)  =  restart  flag. 
Note:  PK(1)  is  equi¬ 

valent  to  Z (151)  and 
in  the  SETUP  deck  must 
be  set  equal  to  PROR 
(Z(l))  .  PK(2)  and 

PK(3)  are  defined  in 
the  restart  deck. 


Used  in  HPHASE  to 
store  the  pressure 
at  the  boundary  be¬ 
tween  cell  k  and  the 
cell  on  its  left,  an 
average  of  the  cell- 
centered  pressures 
of  the  two  cells. 
(Saved  values  of  PRR 
when  column  on  left 
was  being  processed.) 
Storage  used  for 
other  purposes  in 
CDT ,  MAP,  and  TPHASE . 


Flag  set  to  1  when 
using  plugging  failure 
model.  If  PLGOPT=0,  a 
plug  will  not  be  gener¬ 
ated  by  the  code. 


When  PLUGON- 1,  the 
vertical  edge  of  the  plug 
package  has  been  extended 
at  least  one  row  into 
the  target.  (The  plug 
package  now  has  a  non¬ 
zero  volume.)  Used  only 
when  plugging  option  is 
being  used. 


The  total  plastic  work 
done  by  each  material 
package . 


r 


•a 


Vari  ab  Le 
Name 

Location 

Dimension 
of  Array 

Units 

Def init ion 

PLWC 

PLSTC 

(IPLGRT, 
IPLGTP- 
IPLGBT+ 1) 

ergs/g 

An  array  which  stores  a 
specific  plastic  work 
increment  for  each  cell 
in  the  plugging  region 
of  the  target.  Used 
only  when  the  plugging 
option  is  activated  to 
update  the  PLWP  array. 

PLWMIN 

=  Z (59) 

ergs/g 

The  specific  plastic 
work  of  the  material 
near  the  top  right 
edge  of  the  plug  must 
exceed  PLWMIN  before 
the  plug  edge  will  be 
extended . 

PLWP 

PLSTC 

(NTCC) 

ergs/g 

An  array  which  stores  the 
accumulated  specific 
plastic  work  of  the 
material  associated  with 
each  passive  tracer 
particle  in  the  plugging 
region  of  the  target. 

Used  only  when  the 
plugging  option  is 
activated. 

PM  IN 

=  Z (86) 

— 

dynes/cm 

A  pressure  cutoff.  When 
| P(K) | <PMIN  then  P(K)  is 
set  to  zero. 

POSX 

(MOVTCR) 

— 

cm 

The  new  x-coordinate  of 
a  tracer  in  centimeters. 

POSY 

(MOVTCR) 

— 

cm 

The  new  y-coordinate  of 
a  tracer  in  centimeters. 

The  first  two  words  of 
a  restart  dump  are  read 
from  the  file  and  stored 
in  PR(l)(the  flag  word 
=  555.0  or  666.0)  and 


PR 


B.C. 


(IMAX) 


Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


PRCNT 


:Z (16) 


PRDHLT 


:  Z  (45) 


sec 
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of 


PR(2)(the  cycle  number 
the  dump).  In  MAP  the 
PR  array  is  used  to  store 
a  row  of  hollerith  symbols 
to  be  printed  in  one  of 
the  symbolic  maps. 


Convergence  requirement 
for  equilibrating 
pressures  in  a  multi¬ 
material  cell. 


If 


P. 

l 


-  P 


P 


< PRCNT  for 


all  materials  (i) 
cell  K,  then  P(K) 


l  n 


P. 


INPUT  parameter.  Gives 
the  initial  time  inter¬ 
val  between  EDIT  prints. 
There  are  five  parameters 
which  control  printing 
frequency:  PRDELT, 

IPCYCL,  PRLIM,  PREACT, 
and  NUMSCA.  If  the  user 
is  printing  on  time 
(PRDELW  and  IPCYCL=0)  , 
DT  will  be  adjusted  so 
that  a  print  will  occur 
exactly  every  PRDELT 
seconds.  If  the  user  is 
printing  on  cycles 
(PRDELT=0 ,  I PCYCL^O) ,  a 
print  will  occur  every 
IPCYCL  cycles.  PRLIM, 
PR^ACT  and  NUMSCA  are 
used  to  increase  the 
print  interval.  PRLIM 
is  the  time  (or  cycle) 
at  which  PRDELT  (or 
IPCYCL)  and  PRLIM  are 
multiplied  by  PRFACT. 

The  new  value  of  PRLIM 
establishes  the  next 
time  (or  cycle)  when 
the  print  interval  will 
again  be  rescaled.  This 


Variable  Dimension 

Name  Location  of  Array  Units  Definition 


process  continues  at 
n’ost  NUMSCA  times. 

EXAMPLE:  You  wish  to 

print  every  lxlO‘a  sec 
until  you  reach  1x10" 7 
sec,  then  every  1x10' 7 
sec  until  lxlO"6  sec  and 
every  lxlO"6  sec  there¬ 
after: 

PRDELT  =1  x  10' 0 
PRLIM  =1  x  10' 7 
P REACT  =  10 
NUMSCA  =  2. 

PRESUR  B.C.  -  dynes/cm^  Defined  in  EQST :  pressure 

=  f(p,E).  Used  in  CDT  to 
define  P(K)  in  the  case 
of  pure  cells,  and  in  the 
case  of  multimaterial 
cells  to  define  PRS(l,i), 
the  pressure  of  material 

(i). 

PRFACT  =Z(46)  -  -  INPUT  parameter.  Used 

in  EDIT  for  rescaling 
PRDELT,  IPCYCL  and 
PRLIM  when  the  PRLIM 
time  (or  cycle)  is 
reached  (see  PRDELT) . 

Must  be  >  1. 

PRLIM  =  Z  (44)  -  -  INPUT  parameter:  time 

or  cycle  at  which  to 
rescale  PRDELT  (or 
IPCYCL)  and  PRLIM  by 
PRFACT  (see  PRDELT). 

PROB  =  Z(1)  -  -  INPUT  parameter. 

Identifying  problem 
number.  Must  be  be¬ 
tween  0.0001  and 
99.9999,  inclusively. 


7'  x3h*k  vrvr  ->  wjrwjTOSWOTyr 

>2 


Vari ab le 
Name 


Location 


l)i  mens  ion 
of  Array 


Units 


L)c  f  ini  t  i  on 


(HPHASE) 


dynes/cni  The  pressure  at  the 

boundary  between  cell  k 
and  the  cell  on  its 
right,  an  average  of 
the  cell-centered  pres¬ 
sures  of  the  two  cells. 


PRTIME 


=  Z (131) 


PSUM 


(CUT) 


PVOL 


PLSTC 


(IPLGRT) 


PVOL  3 


(PLGMAS) 


PVRTEX 


dynes/cm' 


Initially  set  to  PRDELT 
m  INPUT.  Thereafter 
calculated  in  EDIT. 

When  T= PRTIME,  it  is 
time  to  print  and 
PRTIME  becomes  T  + 

PRDELT,  the  next  time 
to  print . 

Used  as  a  flag  in  CDT. 

If  the  pressures  of  all 
the  materials  in  a  multi¬ 
material  cell  are  less 
than  PMIN ,  PSUM-0,  which 
flags  the  code  to  set 
the  cell  pressure  to 
zero. 

The  volume  asso* 
ciated  with  the  plug 
package  in  a  cell  into 
which  the  vertical  edge 
of  the  plug  has  been  ex¬ 
tended.  Used  to  convert 
target  mass  into  plug 
mass  in  PLGMAS.  Com¬ 
puted  in  PLGVOL.  Used 
only  when  the  plugging 
option  is  activated. 

The  volume  of 
the  target  material  in 
a  cell  into  which  the 
plug  edge  has  been 
extended. 

The  hydrostatic  pressure 
of  the  cell  which  con¬ 
tains  the  vertex  point 
of  the  void  closing  re¬ 
gion.  See  Section  8.4. 
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?S! 

Variable 

Name 


Dimension 

Location  of  Array 


Units 


Def ini tion 


(HPHASE) 


dynes/ cm' 


RATIO 


(CDT, 
DMADJ , 
HPHASE) 


sec 

(in  CDT) 


RELERR 


(EDIT, 

ENCHCK) 


=  Z  (95) 


MXCELL 


(NVOID, 

NMXCLS) 


g/cnT 


The  linear  artificial 
viscosity  term  added  to 
the  top,  bottom,  left, 

and  right  cell  boundary 
pressures,  respectively. 
A  different  viscosity 
can  be  used  by  redefin¬ 
ing  these  variables. 

(See  Section  2. 2. 2. 5) 

Used  in  the  calcula¬ 
tion  of  DT:  the  ratio 

of  (DX,DY)  .  to 
’  'mm 


[(U,V)  + 

L  ^  *  'max 

sound  speed] 


+  local 
dl  _  .  for  a 


- *'~-~Jmax  “ 

given  cell.  In  DMADJ 
and  HPHASE  used  as 
ratio  of  masses  and 
compressions,  respectively 

Used  for  storing  and 
printing  maximum 
relative  error  in  the 
energy  sum. 


A  flag  used  to  activate 
a  grid  rezone.  Can  be 
set  by  the  user  on  a 
restart  cycle  or  by 
the  code  if  NUMREZ  is 
greater  than  0.  (See 
Section  8.1.3) 

The  density  of  the 
materials  in  an  in¬ 
terface  cell.  If  RHO 
(NVOID, m)=l  the  inter¬ 
face  cell  K  (where 
m  =  MFLAG(K)  -100)  is 
a  free  surface  cell. 

If  RHO(n,m)>0 ,  the 
interface  of  package 
n  crosses  the  cell, 
and  if  RH0(n,m)  =  0, 
the  n  interface  does 
not  cross  the  cell, 
and  any  material  n  in 
the  cell  must  be 
"evacuated".  The  RHO 
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Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


RHOC 

(SPHASE) 

(  J.MAX+1 , 3) 

RHOIN 

B.C. 

(NMAT) 

RHOMAX 

(CDT) 

(30) 

RHOW  B.C. 


RHOZ  MXCELL  (30) 


RHOZR  (EQST) 


array  is  used  as  a 
flag  as  well  as  a 
material  property. 

(See  Section  4.2) 

3 

g/cm  The  average  cell 

density  of  three 
rows  of  cells. 

3 

g/cm  The  initial  material 

density  for  each 
package.  Defined  by 
input  cards  read  by 
SETUP. 

3 

g/cm  The  density  of  a 

high  explosive  cor¬ 
responding  to  a 
maximum  point  in 
the  JWL  equation  of 
state.  The  density 
of  HE  in  a  multi¬ 
material  cell  is  not 
allowed  to  exceed 
RHOMAX  during  the 
pressure  iteration. 
Defined  in  a  DATA 
statement  for  four 
explosives . 

3 

g/cm  Density  of  material. 

Defined  in  CDT  and 
used  in  EQST  to 
define  pressure: 

P  =  f (ENERGY, RHOW) . 

g/cm^  Defined  in  DATA 

statement  in  COMDIM, 
the  "included"  element. 
Normal  density  for 
19  materials. 

g/cm  p  of  a  given 

material.  Used  in 
the  Tillotson  and 
JWL  equations  of 
state.  RHOZR=RHOZ (n) 
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Variable  Dimension 

Name  Location  of  Array  Units  Definition 


RMOM  B . C . 


RMU  B.C. 


where  RHOZ  is  de¬ 
fined  in  a  DATA 
statement  in  the 
"included"  element. 
(See  Section  2.3.1.) 

-  g-cm/sec  The  total  post-HPHASE 

radial  momentum  of  a 
slipline  cell.  Com¬ 
puted  in  HPHASE ;  used 
in  UVCALC  to  update 
the  master  and  slave 
material  velocities. 
(See  Section  5.2.) 

(NMAT)  dynes/cm^  Rigidity  modulus  of 

each  material  package. 
Defined  by  input  cards 
read  by  SETUP.  Used 
in  SPHASE. 


ROEPS 


=Z(110) 


INPUT  parameter.  A 
round-off  epsilon  used 
in  calculating  cutoffs 
of  energy,  velocity 
and  mass. 


RTM 


=  Z(57) 


Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
mass  lost  out  right 
side  of  grid. 


RTMU 


=Z(10) 


g-cm/sec 


Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
radial-momentum  lost 
out  right  side  of  grid. 


RTMV 


=  Z (58) 


g-cm/sec 


Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
axial-momentum  lost  out 
right  side  of  grid. 


RTX 


(FRACS) 


The  grid  line  which  is 
the  right  boundary  of 
the  cell  being  processed 
by  FRACS. 
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Variable 

Name 


Dimension 

Location  of  Array 


Uni  ts 


B.C. 


SALPHA 


=  Z ( 1 38) 


radians 


SAMMP 


MXCELL 


(NMAT , 
NMXCLS) 


SAMMU 


(NMAT) 


SAMMV 


(NMAT) 


De  f  i  n  i  t  i  on 

The  distance  from  the 
bottom  of  the  grid  to 
the  intercept  of  an 
interface  with  a  cell's 
right  boundary,  con¬ 
verted  to  cell  units. 

Used  in  NEWR1I0  as  well 

as  in  FRAGS.  (Lquivalenced 
to  WSY  in  FRAGS.) 

The  angle  with  the  x- 
axis  at  which  the  verti¬ 
cal  edge  of  the  plug 
was  last  extended. 
Initially,  SALPHA=tt/2 
in  INPUT.  Used  only 
for  plugging  calcula¬ 
tions. 

The  mass  of  material 
package  n  which  is  trans¬ 
ported  across  the  right 
boundary  of  interface 
cell  K  is  stored  in 
SAMMP (n,m),  where  m  = 

MFLAG (K) - 100 .  If 
SAMMP (n,m)  >  0,  the 
mass  is  moving  out  of 
the  cell;  otherwise, 
it  is  moving  into  the 
cell.  Computed  in 
DMCALC.  (See  Section 
10.4.2) 


g-cm/sec  The  radial  momentum 
of  a  given  material 
transported  across 
the  bottom  of  an 
interface  cell.  (See 
Section  10.4.2) 

g-cm/sec  The  axial  momentum  of 
a  given  material  trans¬ 
ported  across  the 
bottom  of  an  interface 
cell.  (See  Section 
10.4.2) 
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Variable  Dimension 

Name  Location  of  Array  Units  Definition 

SAMMY  MXCELL  (NMAT,  g  The  mass  of  material 

NMXCLS)  package  n  which  is 

transported  across 
the  bottom  boundary 
of  interface  cell  K 
is  stored  in  SAMMY 
(n,m) ,  where  m  = 
MFLAG(K) -100.  If 
SAMMY (n,m)  >0,  the 
mass  is  moving  into 
the  cell;  otherwise, 
it  is  moving  out  of 
the  cell.  Computed 
in  DMCALC.  (See 
Section  10.4.2) 

SAMPY  MXCELL  (NMAT,  .  g  The  mass  of  material 

NMXCLS)  package  n  which  is 

transported  across 
the  top  boundary  of 
interface  cell  K  is 
stored  in  SAMPY (n,m), 
where  m  =  MFLAG(K) 
-100.  If  SAMPY (n,m) 

>  0,  the  mass  is 
moving  out  of  the 
cell;  otherwise ,  it 
is  moving  into  the 
cell.  Computed  in 
DMCALC.  (See  Section 
10.4.2) 

SAMUR  SL  (NMAT)  g -cm/sec  The  radial  momentum 

of  a  given  material 
transported  across 
the  right  of  an  in¬ 
terface  cell.  (See 
Section  10.4.2) 

SAMUT  SL  (NMAT)  g- cm/ sec  The  radial  momentum 

of  a  given  material 
transported  across 
the  top  of  an  inter¬ 
face  cell.  (See 


Variable 

Dimension 

Name 

Location 

of  Array 

Units 

Definition 

SAMVR 

SL 

(NMAT) 

g-cm/sec 

The  axial  momentum 

of  a  given  material 
transported  across 
the  right  of  an 
interface  cell . 

(See  Section  10.4.2) 


SAMVT  SL 


SCALE  (MAP) 


SDELEB  B.C. 


(NMAT) 

g-cm/sec 

The  axial  momentum 
of  a  given  material 
transported  across 

the  top  of  an  in¬ 
terface  cell.  (See 
Section  10.4.2) 

A  scale  factor  used 
to  determine  which 
alphabetic  symbol 
to  associate  with  a 
given  cell  quantity. 

SCALE  is  based  on 
the  minimum  and 
maximum  of  that 
quantity  as  well  as 
the  number  of 
characters  to  be 
used  by  the  symbolic 
map . 

(NMAT) 

ergs 

The  internal  energy 
of  a  given  material 
transported  across 
the  bottom  boundary 
of  an  interface  cell 
(See  Section  10.4.2) 

SDELER 


B.C. 


(NMAT) 


ergs 


SDELET 


B.C. 


(NMAT) 


ergs 
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The  internal  energy 
of  a  given  material 
transported  across 
the  right  boundary 
of  an  interface  cell 
(See  Section  10.4.2) 


The  internal  energy 
of  a  given  material 
transported  across 
the  top  boundary  of 
interface  cell.  (See 
Section  10.4.2) 


Variable 

Name 


Location 


Dimens  ion 
of  Array 


Units 


Definition 


SDELM  B.C. 


SDT  B.C. 


SFLEFT  B.C. 


SGAMC  MXCELL 


SIE  MXCELL 


(NMAT)  g  The  total  change 

in  the  mass  of  a 
given  material  in 
an  interface  cell 
due  to  transport. 

-  sec  Defined  in  INFACE 

(SDT=DT/CYCMX) . 

Time  step  for 
each  subcycle  of 
INFACE. 

(NMAT,  g- cm/sec  The  radial  momentum 

JMAX)  of  a  given  material 

transported  across 
the  left  boundary 
of  an  interface  cell. 
(See  Section  10.4.2) 

(NMAT,  g  The  mass  of  material 

NMXCLS)  package  n  which  is 

transported  across 
the  left  boundary  of 
interface  cell  K  is 
stored  in  SGAMC(n,m), 
where  m  =  MFLAG(K) 
-100.  If  SGAMC (n,m) 

>  0,  the  mass  is 
moving  into  cell  K; 
otherwise,  it  is 
moving  out  of  the 
cell.  Computed  in 
DMCALC .  (See  Section 
10.4.2) 

(NMAT,  ergs/g  The  specific  internal 

NMXCLS)  energy  of  materials 

in  interface  cells. 
These  materials  have 
specific  internal 
energies  different 
from  the  cell  specific 
internal  energy,  AIX, 
which  is  a  mass 
weighted  average  of 
the  material  specific 
internal  energies. 

(See  Section  4.2) 

i  rt  >*(  n  ' 


Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


SIEMIN 

=  Z (82) 

SIGC 

B.C. 

( JMAX) 

SIGMU 

(TPHASE) 

— 

SIGMV 

(TPHASE) 

.... 

SLOPS 

(PLGTCR) 

— 

SLOPT 

(PLGTCR) 

SLPNDX  =Z ( 2  2) 


ergs/g  The  cutoff  on  the 
specific  internal 
energy  increment 
in  the  transport 
phase  (TPHASE) . 

If  this  increment 
is  less  than  SIEMIN, 
the  increment  is 
set  to  zero  and 
ETH  is  adjusted 
accordingly . 

ergs  The  total  energy 

of  the  mass  trans¬ 
ported  across  the 
left  side  of  a  cell 
in  TPHASE.  Equiva- 
lenced  to  PL.  (See 
Section  10.4.2) 

g-cm/sec  Total  change  in  radial 
momentum  of  a  cell. 


g-cm/sec  Total  change  in  axial 
momentum  of  a  cell. 


The  slope  of  the  ex¬ 
tended  edge  of  the  plug. 
See  SLOPT. 

The  slope  of  the  back 
surface  of  the  target 
between  two  consecutive 
tracer  particles.  Used 
to  compute  the  intercept 
of  the  plug  edge  with 
the  back  surface  of  the 
target . 


M 


The  x-coordinate  (in  * 

cell  units)  of  the  top  ? 

of  the  s 1 ipl ine  in  a  c 

plugging  calculation.  j 

Until  the  plug  is  com¬ 
pletely  formed,  the  top 
of  the  si  ip  line  is  the 
last  grid  line  the  fl 
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ar iab le 
Name 

Location 

Dimens  ion 
of  Array 

Units 

Definition 

vertical  edge  of  the 
plug  was  extended  to, 
although  the  top  of  the 
plug  itself  can  move 
beyond  that  grid  line. 
See  TXCL  . 

SLPNDY 

=  Z (23) 

The  y-coordinate  (in 
cell  units)  of  the  top 
of  the  slipline  in  a 
plugging  calculation. 

See  SLPNDX  and  TYCL. 

SNB 

(SPHASE) 

(IMAX+1) 

dynes/cm1" 

Normal  deviator  stress 
at  bottom  of  a  cell . 

SNL 

(SPHASE) 

— 

dynes/cm2 

Normal  deviator  stress 
at  left  of  cell. 

SNLX 

(SPHASE) 

-  -  -  - 

dynes/cm 

=  SNL*X  ( I  - 1 )  for  cell 
in  column  I. 

SNR 

(SPHASE) 

— 

dynes /cm 2 

Normal  deviator  stress 
at  right  of  a  cell. 

SNT 

(SPHASE) 

— 

dynes /cm 2 

Normal  deviator  stress 
at  top  of  a  cell. 

SOLID 

(STRNG , 
SPHASE) 

g/cm3 

=  RHOZi*AMDMi.  If 

RHCK  <. SOLID  for  any 

material  in  cell  K,  the 
deviator  stresses  of 
cell  K  are  set  to  zero. 

SRATIO 

(CDT) 

— 

sec 

Used  to  calculate  DT. 

(SPHASE) 


(IMAX+1,3)  dynes/cm^ 


The  smallest  ratio 
in  the  grid  of  a  cell's 
minimum  dimension  to 
the  sum  of  its  maximum 
velocity  and  sound 
speed. 

The  normal  radial  stress 
deviator  associated 
with  cell  centers.  (See 
Section  2, 3. 2. 2) 
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Variable 

Name 


Dimension 

Location  of  Array  Units 


Definition 


2 

SRRINT  (SPHASL)  -  dynes/cm  The  normal  radial  stress 

deviator  of  the  material 
that  will  be  at  the  cen¬ 
ter  of  cell  K  at  the  end 
of  the  time  step. 

SRZ  (SPHASE)  ( I MAX  +  1,3)  dynes/cm2  The  shear  stress  asso¬ 

ciated  with  cell  centers. 
(See  Section  2. 3. 2. 2) 

The  shear  stress  of  the 
material  that  will  be 
at  the  center  of  cell K 
at  the  end  of  the  time 
s  tep . 

Temporary  storage  for 
the  material  specific 
internal  energies  of 
the  two  cells  being 
combined  into  one  in 
COMPRS. 


SSI  EN 

B.C. 

(NMAT ) 

ergs/g 

Initial  specific  inter¬ 
nal  energy  of  each  material 
package.  Defined  by 
input  cards  read  in 

SETUP . 

SSIGC 

B.C. 

(NMAT, 

JMAX) 

ergs 

The  internal  energy  of 
a  given  material  trans¬ 
ported  across  the  left 
boundary  of  an  inter¬ 
face  cell. 

SSIGMU 

SL 

(NMAT) 

g  -cin/sec 

The  total  change  in  the 
radial  momentum  of  a 
given  material  in  an 
interface  cell  due  to 
transport . 

SSIGMV 

SL 

(NMAT) 

g -cm/sec 

The  total  change  in  the 
axial  momentum  of  a 
given  material  in  an 
interface  cell  due  to 
transport . 

SRZ 1 NT  (SPHASE) 


dynes /cm' 


SSIE  B.C. 


(2,NMAT)  ergs/g 
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Variable  Dimension 

Name  Location  of  Array  Units  Definition 


STAB  =  Z (139) 


STB 

(SPHASE) 

STEZ 

B.C. 

STK1  B.C. 


STK2  B.C. 


STL  (SPHASE) 


(IMAX+1)  iynes/cm^ 

(NMAT)  ergs/g 

(NMAT)  dynes/cm^ 

(NMAT)  dynes/cm^ 

-  dynes/cm^ 


INPUT  parameter. 

Used  in  CDT.  Initial 
value  of  the 
"stability  fraction" 
for  the  calculation 
of  DT.  If  FINAL  = 

0,  STAB  is  constant. 
Otherwise,  its  value 
changes  from  STAB  to 
FINAL  in  a  geometric 
progression.  [Note: 

DT  =  STAB*SRATIO . ] 

Shear  stress  at 
bottom  of  cell. 

E  for  each  material 
m 

package.  Used  in 
yield-strength  cal¬ 
culation  in  STRNG. 
Defined  by  input 
cards  read  by  SETUP. 
(See  Section  2. 3. 2. 3) 

for  each  material 

package.  Used  in 
yield-strength  cal¬ 
culation  in  STRNG. 
Defined  by  input 
cards  read  by  SETUP. 
(See  Section  2. 3. 2. 3) 

Y2  for  each  material 

package.  Used  in 
yield-strength  cal¬ 
culation  in  STRNG. 
Defined  by  input 
cards  read  by  SETUP. 
(See  Section  2. 3. 2. 3) 

Shear  stress  at  left 
of  cell. 


Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


STLX 

(SPHASE) 

.... 

dynes/ cm 

=  STL*X ( I  - 1 )  for 
cell  in  column  I. 

STR 

(SPHASE) 

— 

dynes/cm2 

Shear  stress  at 
right  of  cell. 

STRENG 

(SPHASE, 

STRNG) 

— 

dynes/cm2 

The  shear  yield 
strength  of  cell 

‘  E 

Y=(Y0+Y,M+Y2y2)Cl-Ei) 

m 

It  is  calculated  in 
subroutine  STRNG  and 
stored  in  WS,  a 
blank  common  variable 
then  WS  is  used  to 
define  STRENG  in 
SPHASE. 


STRSRR 

ELPL 

(KMAX) 

dynes/cin  The  cell  centered 
normal  deviatoric 
stress  in  the  radia 
direction . 

STRSRZ 

ELPL 

(KMAX) 

2 

dynes/cm  The  cell  centered 
deviatoric  shear 
stress . 

STRSZZ 

ELPL 

(KMAX) 

2 

dynes/cm  The  cell  centered 
normal  deviatoric 
stress  in  the  axial 
direction . 

STT 

(SPHASE) 

— 

2 

dynes/cm  Shear  stress  at  top 
of  cell. 

SUME 

(TPHASE) 

— 

ergs  Used  in  TPHASE  to 

sum  energy  increments 
being  ignored  on  a 
given  cycle.  (See 
SIEMIN) .  Used  to 
adjust  ETH,  the 
theoretical  energy 
total . 
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Variable 

Name 


Dimension 

Location  of  Array 


Units 


Definition 


TFLUX 


(DMADJ) 


THET 


B.C. 


(2) 


radians 


THETA 


MXCELL 


(NMXCLSJ 


radians 


TH03 


TKO 


(SPHASE) 

(SPHASE) 


1/sec 


dynes2/cm^ 


TK1 


(SPHASE) 


,  2,4 

dynes  /cm 
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The  sum  of  the  mass 
transport  terms  for 
which  the  mass  is  mov¬ 
ing  out  of  a  cell  and 
is  moving  into  a  cell’ 
which  is  not  flagged 
for  evacuation.  These 
"legitimate"  transport 
terms  are  used  to  evacu¬ 
ate  a  material  when  its 
interface  leaves  a  cell. 


Temporary  storage  for 
the  angle  of  the  slip 
line  across  the  two 
cells  being  combined 
into  one  in  COMPRS. 


The  angle  of  the  slip- 
line  (with  reference 
to  the  x-axis)  as  it 
crosses  a  slipline 
cell . 


■  1/3  (ux  +  vy  * 


(See  Section  2. 3. 2.1) 


r2 


=  2Y“,  where 


Y  is  the 
shear  yield  strength 
of  the  material  . 


Twice  the  second 
invariant  of  the  stress 
deviator  tensor: 


s2  +S2  +S2  +2S2  where 
zz  rr  00  rz 


00 


=  -(S  +S  ) 
v  rr  zz 


is  the  hoop  stress. 
Tested  against  TKO 
to  see  if  material 
has  yielded. 


1  I 


Variable 

Name 

Location 

Dimension 
of  Array 

Units 

Definition 

TOPM 

=  Z  (63) 

g 

Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
mass  transported  across 
top  of  grid. 

TOPMU 

=  Z(9)  . 

g- cm/sec 

Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
radial -momentum  trans¬ 
ported  across  top  of 
grid. 

TOPMV 

=Z (66) 

g- cm/sec 

Calculated  in  TPHASE. 
Printed  in  EDIT.  Total 
axial-momentum  trans¬ 
ported  across  top  of 
grid. 

TPX 

B.C. 

— 

- - 

The  distance  from  the 

left  boundary  of  the 
grid  to  the  intercept 
of  an  interface  with 


I  a  cell's  top  boundary. 

Converted  to  cell 
units.  Used  in  NEWRHO 
as  well  as  FRACS. 

(Equivalenced  to  WSX 
in  FRACS.) 

TPY  (FRACS)  -  -  The  grid  line  which 

is  the  top  boundary 
of  the  cell  being 
processed  by  FRACS. 


TRIAL 


(CDT) 


TSTOP 


=Z (50) 


cm/sec  Used  in  CDT.  Maximum 
in  the  grid  of  a 
cell's  sound-speed 
plus  its  largest 
velocity  component. 
Printed  in  CDT  as 
MAX CUV. 


sec 


INPUT  parameter. 
Value  of  T  at  which 
execution  stops  when 
stopping  on  time 
rather  than  cycles. 


Variable  Dimension 

Name  Location  of  Array  Units 


Definition 


TWOPI 

B.C. 

— 

— 

=  2  7T  . 

TX 

B.C. 

(NVOID, 

NTPMX) 

X-coordinates  of  tracer 
particles  circumscribing 
the  material  packages. 

TXCL 

B.C. 

The  x-coordinate  (in 
cell  units)  of  the  top 
right  corner  of  the 
plug  package. 

TX1 

(FRACS) 

The  x-coordinate  (in 
cell  units)  of  the  first 
of  the  pair  of  tracers 
being  considered  by 

FRACS. 

TX2 

(FRACS) 

The  x-coordinate  (in 
cell  units)  of  the 
second  of  the  pair  of 
tracers  being  considered 
by  FRACS. 

TY  . 

B.C. 

(NVOID, 

NTPMX) 

- - 

Y-coordinates  of  tracer 
particles  circumscribing 
the  material  packages. 

TYCL 

B.C. 

The  y-coordinate  (in 
cell  units)  of  the  top 
right  corner  of  the 
plug  package. 

TY1 

(FRACS) 

The  y-coordinate  (in 
cell  units)  of  the  first 
of  the  pair  of  tracers 
being  considered  by 

FRACS. 

TY2 

(FRACS) 

The  y-coordinate  (in 
cell  units)  of  the 
second  of  the  pair  of 
tracers  being  considered 
by  FRACS. 

IJ 

B.C. 

(KMAX) 

cm/sec 

Radial  velocity  of  cell. 

UAMMP 

(TPHASE) 

— 

cm/sec 

Radial  velocity  of  mass 
transported  across 
right  cell  boundary. 
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UAMPY 

(TPHASE) 

- -  cm/ sec 

Radial  velocity  of  mass 
transported  across  top 
cell  boundary. 

(JEFF 

MOV 

-  cm/ sec 

The  area-density  weighted 
average  of  cell  centered 
radial  velocities  used 
to  move  a  given  tracer 
particle.  If  a  tracer 
is  on  the  plug  package 
boundary,  UEFF  is  com¬ 
puted  in  PLUGUV. 

UK 

(S PHASE) 

(IMAX+1 , 3)  cm/sec 

Temporary  storage  for 
three  rows  of  radial 
velocities . 

UKT 

(S PHASE) 

-  cm/sec 

Temporary  storage  for 

U(K)  when  computing  the 
elastic-plastic  work 
done  by  cell  K. 

UL 

B.C. 

(See  definition  cm/sec 
of  JDX2  param¬ 
eter  in 

Appendix  A) 

* 

----  cm/sec 

Radial  velocity  at  left 
boundary  of  cell  K, 

Used  in  HPHA5E.  Storage 
for  other  purposes  in 

CDT,  MAP,  and  TPHASE. 

UMIN 

=  Z  (129) 

Calculated  in  CDT. 

Used  as  velocity  cutoff 
in  TPHASE,  DMCALC ,  and 
MOVTCR. 

UMS 

(UVCALC, 

UVMOD) 

-  cm/sec 

A  mass  weighted  average 
of  the  radial  velocities 
of  all  the  "master” 
materials  in  a  given 
slipline  cell. 

UMSNEW 

(UVCALC) 

cm/ sec 

The  post-HPHASE 
radial  velocity  of  the 
"master"  materials  in 
a  slipline  cell.  (See 
Section  5.2) 

UNxxx 

=Z (xxx) 

-  - 

Unused  Z-storage. 

10-83 


m 


'•1 


Variab le 
Name 


Location 


Dimension 
of  Array 


Units 


Definition 


URR 

R.C. 

cm/sec 

The  radial  velocity  at 
the  right  boundary  of 
cell  K,  an  average  of 
the  cell-centered  velo¬ 
cities  of  the  two  cells. 
Used  in  HPIIASH,  DMCALC 
and  TPHASE . 

US 

MXCELL 

(NMAT , 
NMXCLS) 

cm/sec 

The  radial  velocity  of 
a  material  in  an  inter¬ 
face  cell.  The  mater¬ 
ials  in  a  slipline  cell 
have  different  radial 
velocities.  The  cell- 
centered  radial  velocity 
(U)  of  a  slipline  cell 
is  a  mass -weighted 
average  of  the  material 
radial  velocities  (US). 
(See  Section  4.2) 

USL 

(UVCALC, 

IJVMOU) 

cm/sec 

A  mass  weighted  average 
of  the  radial  velocities 
of  all  the  "slave" 
materials  in  a  given 
slipline  cell. 

USLNEW 

(UVCALC) 

/ 

cm/sec 

The  post-HPHASE  radial 
velocity  of  the  "slave" 
materials  in  a  slipline 
cell.  See  Section  5.2. 

USSAVE 

B.C. 

(2,WT) 

cm/sec 

Temporary  storage  for 
the  material  radial 
velocities  of  two  cells 
being  combined  into  one 
in  COMPRS. 

UTK 

MOV 

(4) 

cm/ sec 

The  radial  velocity 
components  of  the 
four  overlap  cells  used 
to  compute  the  radial 
velocity  of  a  tracer 
particle . 

UTL 

B.C. 

(NMAT , JMAX) 

cm/sec 

The  radial  velocity 
of  a  given  material 
being  transported 
across  the  left  boundary 

of  a  cell. 
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Variable  Dimension 

Name  Location  of  Array  Units  Definition 


U TRANS 

B.C. 

(NMAT , 4) 

cm/ sec 

The  radial  velocity  of 
a  given  material  trans¬ 
ported  across  one  of 
the  four  boundaries  of 
cell  K  (l=top,  2=bottom, 
3=right ,  4=left) . 

UUR 

B.C. 

(NMAT) 

cm/sec 

Initial  radial  velocity 
of  each  material  pack¬ 
age.  Defined  by  in¬ 
put  cards  read  in  SETUP. 

V 

B.C. 

(KMAX ) 

cm/sec 

Axial  velocity  of  cell. 

VABOVE 

B.C. 

cm/sec 

The  axial  velocity  at 
the  top  boundary  of 
cell  K ,  an  average  of 
the  cell  centered 
velocities  of  the  two 
cells.  Used  in  HPHASE , 
DMCALC  and  TPHASE. 

VALUE 

(MAP) 

(41) 

-  -  -  - 

Maximum  value  associated 
with  each  symbol  print¬ 
ed  by  MAP. 

VAMMP 

(TPHASE) 

cm/sec 

Axial  velocity  of  mass 
transported  across  right 
cell  boundary. 

VAMPY 

(TPHASE) 

' 

cm/sec 

Axial  velocity  of  mass 
moving  across  top  cell 
boundary. 

VBLO 

(H PHASE) 

— 

cm/sec 

The  axial  velocity  at 

the  bottom  boundary  of 
cell  k,  an  average  of 


the  cell-centered  velo¬ 
cities  of  the  two  cells. 

(Local  -  cm^  The  total  volume  of  a 

variable  cell.  VCELL-TAU(I)*DY(J) . 

in 

several 

subroutines) 
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VCELL 


Variable 

Name 


Dimension 

Location  of  Array 


Un  i  ts 


Def inition 


VCT 

(DMCALC) 

0 

cm 

The  minimum  of  the 
volume  of  the  donor 
and  acceptor  cells. 

Calculated  and  used  in 

DMCALC  to  test  if  a 
mass  transport  term  is 
too  small  to  be  counted. 

VEFF 

MOV 

cm/sec 

The  area-density  weighted 
average  of  ce 11 -centered 
axial  velocities  used  to 
move  a  given  tracer  par¬ 
ticle.  If  a  tracer  is 
on  the  plug  package 
boundary,  VHFF  is  com¬ 
puted  in  PLUGUV. 

VK 

(SPHASE) 

(IMAX+1,3) 

cm/sec 

Temporary  storage  for 
three  rows  of  axial 
velocities . 

VKT 

(SPHASE) 

cm/ sec 

Temporary  storage  for 

V(K)  when  computing  the 

elastic-plastic  work 

done  by  cell  K.  j 

VMS 

(UVCALC, 

UVMOD) 

cm/sec 

A  mass  weighted  average  j 

of  the  axial  velocities 
of  all  the  "master" 
materials  in  a  given 
slipline  cell.  j 

VMSNEW 

(UVCALC) 

cm/ sec 

The  post-UPHASE  axial 
velocity  of  the  "master" 
materials  in  a  slipline  * 

cell.  See  Section  5.2. 

VOLM 

MXCELL 

(NVOID, 

NMXCLS) 

3 

cm 

The  partial  cell  volume 
of  a  material  in  an  ; 

interface  cell.  Com¬ 
puted  by  VOLFND  to  '  : 

determine  the  mass  of  * 

each  material  of  an  in¬ 
terface  cell.  Used 
only  when  generating  a 
problem.  liqui  valenced 
to  SAMMY  array.  (Note: 
this  array  is  larger 
than  the  SAMMY  array 
and  overflows  into  the 
SGAMC  array.) 
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Variable 

Dimens  ion 

N  ame 

Locat ion 

of  Array 

Uni  ts 

Definition 

VOW 

(EQST) 

-  -  -  - 

G. 

o 

a 

ii 

VS 

MXCELL 

(NMAT , 
NMXGLS) 

cm/ sec 

The  axial  velocity  of 
a  material  in  an  inter 

face  cell.  The  mater¬ 
ials  in  a  slipline  cell 
have  different  axial 


velocities.  The  cell- 
centered  axial  velocity 
(V)  of  a  slipline  cell 
is  a  mass-weighted 
average  of  the  material 
axial  velocities  (VS) . 

(See  Section  4.2.) 


VSL 

(UVCALC, 

UVMOD) 

cm/sec 

A  mass-weighted  average 
of  the  axial  velocities 
of  all  the  "slave" 
materials  in  a  given 
slipline  cell. 

VSLNEW 

(UVCALC) 

cm/sec 

The  post-HPHASE  axial 
velocity  of  the  "slave" 
materials  in  a  slipline  ;j 

cell.  See  Section  5.2. 

VSSAVE 

B.C. 

(2  ,NMAT) 

cm/sec 

Temporary  storage  for 
the  material  axial 
velocities  of  two  cells 
being  combined  into 
one  in  COMPRS. 

VSUM 

(CDT , 
FILGRD , 
TPHASE) 

- - 

3 

cm 

The  sum  of  the  partial 
volumes  of  the  materials 
in  an  interface  cell. 

VTK 

MOV 

(4) 

cm/sec 

The  axial  velocity 
component  of  the  four 
overlap  cells  used  to 
compute  the  axial 
velocity  of  a  tracer 
particle. 

VTL 

B.C. 

(NMAT, JMAX) 

cm/ sec 

The  axial  velocity  of 

a  given  material  being 
transported  across 
the  left  boundary  of 
cell  K. 


Variable 

Dimension 

Name 

Location 

of  Array 

Units 

Definition 

VTRANS 


B.C. 


(NMAT,4)  cm/ sec 


(NMAT) 


cm/,sec 


WFLAGF 


»Z(51) 


WFLAGL 


=  Z  ( 5  2 


WSMAX 


WSMIN 


(MAP) 


(MAP) 


The  axial  velocity  of 
a  given  material 
transported  across  one 
of  the  four  boundaries 
of  cell  K  (l=top,  2= 
bottom,  3=right,  4= 
left). 

Initial  axial  velocity 
of  each  material 
package.  Defined  by 
input  cards  read  by 
SFTIJP . 

The  weighting  factor 
associated  with  each  of 
the  four  overlap  cells 
used  to  compute  the 
velocity  components  of 
a  tracer  particle. 

Used  by  MOVTCK  and 
PLUGUV. 

Used  in  INPUT  and  EDIT. 
Set  =  1  on  first  cycle 

of  a  run  in  INPUT. 
Triggers  an  EDIT  print 
on  first  cycle  of  every 
run.  Reset  to  0  at  end 
of  EDIT. 

Used  in  HELP  and  EDIT. 
Flags  last  cycle.  Set 
1  in  EDIT.  Signals 
HELP  to  call  exit. 


Used  in  most  sub¬ 
routines  for  working 
storage . 

The  maximum  value  of 
each  quantity  repre¬ 
sented  by  the  symbolic 
maps . 

The  minimum  value  of 
each  quantity  represent 
ed  by  the  symbolic  maps 
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Variable 

Name 


Location 


Dimension 
of  Array  Units 


Definition 


WSQR  B.C.  (NMAT) 


WSUM  (CDT, 

PLGGEN, 

SPHASE) 


X  B.C.  (IMAX) 


X00  B.C. 


XAXIS  MISC  (NVOID, 

IMAX) 

\ 


XIENRG  =Z (140) 


cm  /  A  term  in  the  weight- 

dynes-g  ing  factor  used  in 
the  iteration  to 
equilibrate  material 
pressures  in  multi¬ 
material  cells: 

WSQRi=l/(3Pi/8Vi)E  = 

1/(P?CSQR.)  . 

(See  Section  4.6.1) 

| 

dynes?  The  sum  of  the  weight- 

g  /cnr  ing  factors  used  in 

the  pressure  iteration: 

WSUM  -J^/WSQR.  . 

(See  Section  4.6.1) 

In  PLGGEN  and  SPHASE 
WSUM  is  the  sum  of  the 
area  weighting  factors. 

cm  The  distance  from  the 

left  grid  boundary  to 
the  right  edge  of  a 
cell. 

cm  The  x-coordinate  of 

X (0) ,  the  left  grid 
boundary  (assumed 
to  be  zero) . 

-  The  x-coordinate  (in 

cell  units)  of  a 
tracer  on  the  x-axis 
of  the  grid.  If 
XAXIS(n,I)>0  the  inter¬ 
face  of  package  n  has 
a  tracer  on  the  x-axis 
in  column  I. 

ergs  Total  internal  energy 

in  the  grid.  Calculated 
in  EDIT  and  used  for 
printing  labels  on 
tracer  particle  plots. 


XKENRG  =Z(141)  -  ergs  Total  kinetic  energy 

in  the  grid.  Calcu¬ 
lated  in  EDIT  and 
used  for  printing 
labels  on  tracer 
particle  plots. 

XMAS  B.C.  (2fNMAT)  g  Temporary  storage 

for  the  material 
masses  of  two  cells 
being  combined  into 
one  in  COMPRS. 


XMASS 

MXCELL 

(NMAT , 
NMXCLS) 

g 

The  mass  of  a 
material  in  an  inter¬ 
face  cell.  (See 
Section  4.2) 

XMMS 

(IJVCALC , 
UVMOD) 

g 

The  sum  of  the  masses 
of  the  "master" 
materials  in  a  slip¬ 
line  cell. 

XMSL 

(UVCALC, 

UVMOD) 

g 

The  sum  of  the  masses 
of  the  "slave" 
materials  in  a  slip¬ 
line  cell. 

XMSM 

(PLGMAS) 

g 

The  sum  of  the  plug 
and  target  masses  in 
a  cell  in  which  all 
the  target  mass  is 
being  converted  to 
plug  mass.  (See 
Section  6.2.4) 

XP 

TRACRS 

(NTCC) 

“  “  “  ” 

X-coordinates  of 
passive  tracer 
particles . 

XRH 

B.C. 

(2, NMAT) 

g/cm3 

Temporary  storage 

for  the  material 
density  of  two  cells 
being  combined  into 
one  in  COMPRS. 
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Variable 

Name 


Location 


Dimension 
of  Array 


Units 


Definition 


XSAVE  (FRACS) 


XTENRG  =  Z(142)  -  ergs 


XUM  (MAP)  (41) 


Y  B.C.  (JMAX)  cm 

Y00  B.C.  -  cm 


The  x-coordinate  (in 
cell  units)  of  the 
previous  intercept 
of  an  interface  with 
a  given  cell  on  the 
x-axis.  Covers  the 
case  where  an  inter¬ 
face  begins  and  ends 
in  the  same  cell  on 
the  x-axis. 

Total  energy  in  the 
grid.  Calculated  in 
EDIT  and  used  for 
printing  labels  on 
tracer  particle  plots. 

Used  in  MAP.  Array 
has  negative  alphabetic 
characters  for  symbolic 
maps.  Defined  in  a 
DATA  statement. 

The  distance  from  the 
bottom  of  the  grid  to 
top  of  a  cell. 

The  y-coordinate  of 
Y(0),  the  bottom  grid 
boundary  (assumed  to 
be  zero) . 


YAMC 


B.C. 


(JMAX) 


g-cm/sec 


Calculated  and  used 
in  TPHASE.  Axial  momentum 
of  mass  transported 
across  left  side  of 
cell.  Equivalenced  to 
UL  array.  (See  Section 
10.4.2) 


YAXIS 


MISC 


(NVOID, 

JMAX) 


The  y-coordinate  (in 
cell  units)  of  a 
tracer  on  the  y-axis 
of  the  grid.  If 
YAXIS (n, J) >0 ,  the 
interface  package  n 
has  a  tracer  on  the 
y-axis  in  column  J. 
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Variab le 
Name 


YSAVE 


Dimension 

Location  of  Array 
( FRACS }  - 


Uni  ts 


ZMOM 


Definition 

The  y-coordinate  (in 
cell  units)  of  the 
previous  intercept  of 
an  interface  with  a 
given  cell  on  the  y- 
axis.  Covers  the  case 
where  an  interface  be¬ 
gins  and  ends  in  the 
same  cell  on  the  y-axis. 

Storage  for  most  of  the 
input  parameters  follows 
Z,  the  first  word  of 
blank  common.  The  "Z- 
block"  (the  first  150 
words  of  blank  common) 
is  written  on  a  file  for 
restarting  problems. 


g-cm/sec  The  total  post  - Hl’llASli 
axial  momentum  of  a 
slipline  cell.  Com¬ 
puted  in  HPUASIi;  used 
in  UVCALC  to  update 
the  master  and  slave 
material  velocities. 
See  section  5.2. 


10.4  FLAGS  AND  CONVENTIONS 


10.4.1  Flags  Governing  Interface  Cells 


FORTRAN  Statements 


Meaning 


RHO (1 »M)  =  -1.0 


NVOID  =  NMAT  +  1 


M  =  MFLAG(K)  -  100 
RHO (NVOID, M)  =1.0 

MFLAG(K)  =  0 


The  M  location  in  the  material  arrays 
(XMASS,  RHO,  SIE ,  FRACTP ,  FRACRT ,  US, 
VS,  THETAS,  SAMMP,  SAMPY,  SAMMY, 
SGAMC)  is  not  in  use. 

NVOID,  the  number  of  the  "void" 
package,  is  always  one  greater  than 
the  number  of  material  package:;. 

Cell  K  contains  the  free  surface 
interface . 

Cell  K  is  empty  and  is  not  cut  by  an 
interface . 


0  <  MFLAG(K)  <  100 


MFLAG(K)  =  2 


MFLAG(K)  >  100 


Cell  K  is  nonempty  and  pure;  i.e.,  it 
does  not  contain  an  interface. 

Cell  K  is  completely  inside  the 
material  package  2  boundary  and 
contains  only  package  2  material. 

Cell  K  is  an  interface  cell;  i.e.,  it 
contains  at  least  one  interface.  M  = 
MFLAG(K)  -  100  gives  location  of 
quantities  in  the  material  arrays  for 
cell  K. 


MFLAG(K)  <  0 


MAT (2)  =  3 


Cell  K  was  an  interface  cell  on  the 
previous  cycle  but  is  no  longer  cut 
by  an  interface  and  will  be  reflagged 
at  the  end  of  TPHASE. 

The  material  code  number  of  package  2 
is  3.  The  list  in  EQST  indicates  that 
code  number  3  corresponds  to  iron. 

The  material  of  package  1  is  an  ideal 
gas  which  is  given  special  treatment  in 
the  sound  speed  calculation  and  in  the 
strength  phase  of  the  code. 


MAT ( 1 )  =  20 


FORTRAN  Statements 


M  =  MFLAG(K)  -  100 
RHO  (N  ,M)  =  0. 


M  =  MFLAG(K)  -  100 
RHO (N ,M) >  0 . 

XMASS (N  ,M)  =  0. 


M  =  MFLAG(K)  -  100 
RHO  (N  ,M)  =  0. 

XMASS (N,M)  >  0. 


M  =  MFLAG(K)  -  100 
THETA  CM)  -  -  1.0 


M  =  MFLAG(K)  -  100 1 
THETA (M)  >0.  ) 


' A '•f’l  C-  xm  .* 


Meaning 


The  interface  of  package  N  does  not 
cut  the  interface  cell  K.  Material  of 
package  N  will  not  be  transported  into 
cell  K. 


The  package  N  interface  cuts  the  inter¬ 
face  cell  K,  but  cell  K  does  not  contain 
any  material  of  package  N.  Either  no 
material  N  has  yet  entered  the  cell  or 
it  has  all  been  transported  out  before 
the  package  N  interface  has  left  the 
cell. 


The  package  N  interface  has  left  cell 
K,  and  the  material  of  package  N 
[ XMASS (N,M)]  should  be  completely 
evacuated  from  cell  K  on  this  cycle. 


A  slipline  does  not  cross  cell  K,  and 
the  material  velocities,  (US,  VS)  are 
the  same  as  the  cell  velocities,  (U,  V)  . 


A  slipline  crosses  cell  K  at  an  angle 
with  the  x-axis  equal  to  THETA(M) 
radians.  The  material  velocities 
(US,  VS)  of  the  master  and  slave 
packages  will  not  (in  general)  be  equal, 
and  the  cell  velocities  (U,  V)  will  be 
a  mass -weighted  average  of  the  material 
velocities.  See  Chapter  V. 
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10.4.2  Definition  of  Transport  Variables 

In  HELP  there  are  four  transported  quantities:  mass, 
radial  momentum,  axial  momentum  and  total  energy  (or 
internal  energy  in  the  case  of  interface  cells) .  There  are 
sixteen  variables  which  store  these  four  quantities  for  the 
four  boundaries  of  a  cell;  these  variables  are  indicated  in 
the  Figure  10.2.  The  direction  of  the  arrows  indicates 
the  direction  of  flow  for  which  the  variables  are  positive. 
Negative  values  of  these  variables  indicate  the  flow  is 
opposite  to  the  direction  of  the  arrows. 


MASS 


AMPY 


AMMY 


RADIAL  MOMENTUM 


AMUR 


AMMU 


ENERGY 


DGEET 


sigc 


DEI, HR 


DEL  Eli 


AXIAL  MOMEN' [ UM 
AMY  r 


YAMC 


AMY'R 


AMMV 


Figure  10.2--The  sixteen  variables  which  store  the 

transported  quantities  between  pure  cells, 
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The  variables  given  in  Figure  10.2  are  cell  variables,  and 
are  used  to  compute  the  effects  of  transport  for  a  pure  coll. 

An  interface  cell,  on  the  otherhand,  must  have  each  of 
these  transport  quantitites  defined  for  each  material  package. 
The  material  transport  variables  were  named  by  preceding  the 
corresponding  cell  transport  variable  names  with  an  "S".  For 
example,  the  mass  of  package  2  transported  across  the  top  of 
cell  K  is  stored  in  SAMPY  (2,  m) ,  whereas  the  total  mass 
transported  across  the  top  of  cell  K  is  stored  in  AMPY. 
Furthermore,  the  cell  variables  of  interface  cells  are  defined 
by  summing  the  material  variables.  For  example,  for  an  inter¬ 
face  cell 

4MPY  =  2  SAMPY (n,m) 
n 

where  m  =  MFLAG(K)  -  100.  Only  the  interface  cell  energy 
variables  do  not  exactly  correspond  to  the  pure  cell  variables. 
For  example,  DELET  is  the  total  energy  transported  across  the 
top  of  cell  K  and  SDELET  (n,m)  is  the  internal  energy  of 
package  n  transported  across  the  top  of  the  cell.  Therefore, 

DELET  =  £  SDELET(n,m)  +  j  £[ SAMPY (n ,m)  •  (udjj+vd^)  ] 
n  n 

where  ud  and  vd  are  the  donor  cell  material  n  velocity 
n  n  ' 

components.  The  energy  variables  at  the  other  boundaries  are 
defined  in  a  like  manner. 

10.4.3  Radial  and  Axial  Terms 

Both  x  and  r  are  used  interchangeably  in  the  FORTRAN 
variable  names  in  the  code  and  in  the  notation  of  this  report 
to  denote  the  radial  direction.  For  example,  the  radial  di¬ 
mensions  of  the  cells  are  stored  in  the  DX  array,  but  are 
printed  out  under  the  heading  "DR".  The  following  equation 


CHAPTER  XI 


SAMPLE  HELP  INPUT  AND  OUTPUT 

The  HELP  input  and  output  will  be  illustrated  in  this 
chapter  by  considering  three  representative  calculations: 

1.  Shaped  Charge  Liner  Collapse 

2.  Perforation  of  a  Thin  Target 

3.  Impact  into  a  Thick  Target. 

The  input  cards  for  these  problems  will  be  listed, 
and  the  distinctive  aspects  of  each  problem  will  be  discussed. 
The  output  resulting  from  cycle  0  of  each  calculation  will 
also  be  presented,  as  well  as  configuration  plots*  showing  the 
material  interfaces  and  the  calculational  grid.  Some  impor¬ 
tant  aspects  of  the  HELP  output  will  be  discussed  in 
Section  11.4. 

H.].  SHAPED  CHARGE  LINER  COLLAPSE 

The  shaped  charge  calculation  presented  here  will 
illustrate  the  use  of  sliplines  and  high  explosive  detona¬ 
tion.  The  calculation  consists  of  three  material  packages 
and  a  void  (NMAT  =  3):  package  one  is  a  copper  liner; 
package  two  is  a  charge  of  Comp  B,  which  is  detonated  at 
a  point  on  the  axis  of  symmetry;  package  three  is  an 
aluminum  casing;  and  package  four  is  the  void.  (The  void 
package  number  is  always  NMAT  +1.) 

Figure  11.1  shows  the  initial  configuration  of  the 
problem,  and  Figure  11.2  shows  the  initial  zoning  of  the 
grid  relative  to  the  dimensions  of  the  material  packages. 

For  this  problem  the  finest  zoning  was  chosen  to  be  in  the 
region  of  the  grid  where  the  high  explosive  detonates  and 
where  the  jet  forms  since  it  is  important  to  resolve  the 
detonation  front  and  the  stagnation  region  behind  the  jet. 

Because  of  the  axis  of  symmetry  in  these  calculations,  only 
the  right  half  of  each  figure  is  actually  in  the  computational 
grid. 


Figure  11.2 


Initial  zoning  of  the  grid  relative  to  the 
dimensions  of  the  liner,  explosive  and  cas 
in  the  shaped  charge  liner  collapse  calcul 


The  coarest  zoning  is  to  the  right  of  the  casing  since 
this  region  is  not  of  great  interest  in  the  solution. 

Since  this  is  a  shaped  charge  calculation,  the  vari¬ 
able  NLINER  is  set  to  1,  the  liner  package  number.  The 
value  of  the  variable  NOSLIP  remains  0,  its  default  value, 
indicating  that  the  slipline  option  will  be  used.  All  of 
the  interfaces  between  the  explosive  and  the  metal  are 
designated  as  sliplines  since  it  is  assumed  that  the  ex¬ 
panded  detonation  products  do  actually  "slip"  on  the  metal 
surfaces.  Both  the  liner  and  casing  are  defined  to.be 
"master"  packages;  the  explosive  is  defined  to  be  a  "slave" 
package.  (As  discussed  in  Chapter  V,  the  symmetry  of  the 
equations  governing  the  motion  of  "master"  and  "slave" 
packages  makes  the  designation  of  packages  as  "master"  or 
"slave"  an  arbitrary  one.) 

Since  tremendous  distortions  are  expected  for  the 
material  in  the  slug  and  jet  portions  of  the  liner,  the  • 
variables  NADD,  MINX,  MAXX,  MINY,  and  MAXY  are  defined  so 
that  every  ten  cycles  (NADD  =  10)  ADDTCR  checks  tracer 
spacing  in  that  region  and  adds  tracers  when  necessary. 

Because  of  the  geometry  of  the  liner  and  the  high  ex¬ 
plosive,  only  the  portion  of  the  explosive  package  which  is 
shaded  in  Figure  11.3  can  be  "seen"  by  the  primary  initia¬ 
tion  point  (see  Section  7.2.7) .  Thus  it  is  necessary  to 
define  a  secondary  initiation  point,  denoted  by  ®  on  the 
figure,  to  calculate  detonation  times  for  the  unshaded 
part  of  the  explosive. 

If  the  calculation  is  to  be  carried  to  late  times 
after  the  jet  has  formed,  the  user  may  want  to  rezone  the 
grid  in  the  axial  direction  to  keep  the  jet  tip  within  the 
calculational  grid.  This  must  be  decided  at  the  time  the 
problem  is  generated  since  in  order  to  rezone  the  grid, 
JMAX  (the  number  of  rows  in  the  grid)  must  be  an  even 
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number.  Then,  later  in  the  calculation,  just  before  the 
jet  tip  enters  the  top  row  of  the  grid,  the  calculation 
must  be  stopped  and  the  rezone  flags  (REZ  and  JEXTY)  set 
in  the  restart  deck  so  that  the  grid  will  be  rezoned  on 
the  restart  cycle.  The  rezoned  grid  will  have  less 
resolution  in  the  axial  direction,  but  it  will  have 
JMAX/2  more  void  cells  into  which  the  jet  can  expand. 

This  of  course,  is  possible  only  if  the  problem  is 
generated  so  that  the  jet  moves  toward  the  top  rather 
than  the  bottom  grid  boundary.  (See  Section  8.1  for  more 
details  on  7'ezoning  procedures.  In  particular,  see 
Section  8.1.2  for  the  reason  the  jet  tip  cannot  be 
allowed  to  enter  top  row  of  grid  before  the  grid  is 
rezoned . ) 

The  input  cards  for  the  shaped  charge  problem  are 
listed  on  the  following  pages  along  with  the  cycle  0  out- 
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for  Package  1 
(Liner) 
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Material  Tracers 
for  Package  2 
(Explosive) 


Material  Tracers 
for  Package  3 
(Casing) 


Material  Tracers 
for  Package  4 
(Void) 


End  of  Material  Tracer  Data 
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11.2  PERFORATION  OE  A  THIN  TARGET 


This  calculation  is  presented  here  to  illustrate  the 
use  of  the  HELP  plugging  option  to  simulate  the  perforation 
of  a  thin  metal  plate  by  a  blunt  metal  projectile.  (See 
Chapter  VI  on  the  HELP  plugging  model  and  Section  12.2  for 
results  of  a  plugging  calculation.)  As  dictated  by  the 
plugging  model,  this  calculation  consists  of  three  material 
packages  and  a  void  (NMAT  =  3):  package  one  is  the  steel 
projectile  and  its  initial  axial  velocity  is  3.71  x  10^cm/sec 
package  two  is  to  be  the  plug  (the  plug  has  no  volume  at 
t  =  0.);  package  three  is  the  aluminum  target,  and  package 
four  is  the  void.  (The  void  package  number  is  always 
NMAT  +1.) 

Figure  11.4  shows  the  initial  configuration  of  the 
problem  as  well  as  the  passive  tracer  particles  which  are 
automatically  generated  in  the  plugging  region  of  the  target 
when  PLGOPT  -  1.  Figure  11.5  shows  the  initial  zoning  of  the 
grid  relative  to  the  dimensions  of  the  projectile  and  target. 
The  finest  zoning  was  chosen  to  be  in  the  region  of  the 
impact  and  in  the  region  where  the  plug  is  expected  to  form. 
The  value  of  DY(1)  was  chosen  so  that  the  initial  projectile- 
target  interface  would  be  coincident  with  a  grid  line,  a 
desirable  procedure  but  not  a  requirement.  The  coarsest 
zoning  is  at  the  lower  portion  of  the  projectile  and  at  the 
right  end  of  the  target,  where  there  is  no  need  for  a  highly 
resolved  solution. 


11-38 


Figure  11 . 4- - Initial  configuration  of  the  perforation  of  a  thin  target  calculation. 

The  dots  represent  the  passive  tracers  which  are  automatically  generated 
in  the  plugging  region  of  the  target. 
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The  Z-block  variables  which  must  be  defined  for  a 


plugging  calculation  are:' 


Z(lll)  PLGOPT  a  flag  which,  when  set  to  1.0, 

activates  the  plugging  mechanisms 


Z  (59) 


PLWMIN 


the  specific  plastic  work  limit  for 
advancing  the  vertical  edge  of  the 
plug 


Z  ( 56) 
Z  (60) 
Z  (61) 


IPLGRT 


IPLGBT 


IPLGTP 


the  rightmost  column,  the  bottom  row 
and  the  top  row  of  the  region  of  the 
target  where  the  passive  tracers  are 
generated  and  where  the  plug  is 
expected  to  form 


Z(109)  NVRTEX  the  second  index  of  the  free  surface 


tracer  at  the  vertex  of  the  void 


closing  region  (the  top  right  corner 
of  the  projectile) 


Z  (112)  NSLD 


the  maximum  number  of  si ip line  cells 
to  be  generated.  The  plugging  option 
presumes  the  use  of  a  slipline  along 
the  vertical  edge  of  the  plug. 


The  projectile  and  plug  packages  must  be  designated  as 
"slaves"  and  the  target  as  a  "master"  since  the  plugging 
routines  make  these  assumptions.  (Details  on  generating 
the  tracer  particles  and  sliplines  for  a  plugging  calculation 
are  given  in  Section  6.1.)  The  user  should  read  Chapter  VI 
before  attempting  to  apply  the  plugging  option  to  an  impact 
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situation.  In  its  present  form  the  plugging  model  has 
several  limitations  that  the  user  should  be  aware  of  before 
using  it. 


The  input  cards  for  the  plugging  problem  are  listed  on 
the  following  pages  along  with  the  cycle  0  output. 
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11.3  IMPACT  INTO  A  THICK  TARGET 


This  calculation  illustrates  the  impact  of  a  heavy 
metal  projectile  into  a  thick  metal  target.  This  is  a 
situation  in  which  plugging  failure  is  not  expected,  there¬ 
fore  the  plugging  option  is  not  used  and,  by  default, 

PLGOPT  =  0.  In  this  calculation  there  are  two  material 


packages  and  a  void  (NMAT  =  2) :  package  one  is  the  tungsten 

_ ■ '  i  _ 3  ‘  j ■  j  ..  i  „  .  •  „  i  l  j •  _  r  r  ...  i  n  4 / 


projectile  and  its  initial  axial  velocity  is  5.5.x  10  cm/sec; 
package  two  is  the  steel  target;  and  package  three  is  the 
void.  (The  void  package  number  is  always  NMAT  +  1.) 


Figure  11.6  shows  the  initial  configuration  of  the 
problem,  and  Figure  11.7  shows  the  initial  zoning  of  the 
grid  relative  to  the  dimensions  of  the  projectile  and  target. 
The  finest  zoning  was  chosen  to  be  in  the  region  of  the 
impact  where  the  greatest  material  deformation  is  expected 
to  occur.  The  value  of  DY(1)  was  chosen  so  that  the  initial 
projectile-target  interface  would  be  coincident  with  a  grid 
line,  a  desirable  procedure  but  not  a  requirement. 


Sliplines  are  not  necessary  in  this  calculation,  so 
NOSLIP  is  set  to  1,  and  the  cards  defining  the  slipline  end¬ 
points  are  omitted  from  the  input  deck.  The  automatic  void 
closing  routine  is  used,  however.  It  is  activated  by 
setting  NVRTEX  -•  T08  which  indentifies  the  void  tracer  at  the 
vertex  point,  which  is  the  intercept  between  the  target  and 
projectile  edge. 


The  input  cards  for  this  calculation  are  listed  on  the 
following  pages  along  with  the  cycle  0  output. 
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CHAPTER  XII 


REPRESENTATIVE  APPLICATIONS  OF  THE  HELP  CODE 

Over  the  past  several  years  the  HELP  code  has  been 
applied  to  a  wide  variety  of  physical  problems  in  the  areas 
of  fluid  and  solid  mechanics.  These  applications,  which  have 
involved  a  large  variety  of  materials  with  a  wide  range  of 
material  properties  and  initial  conditions,  include  studies 
in  hypervelocity  and  ballistic  impact,  high  explosive-metal 
systems,  shaped  charge  jet  formation  and  penetration,  and 
energy  deposition  and  stress  wave  propagation.  In  this 
section  selected  results  from  these  various  applications  are 
presented.  References  are  provided  for  those  interested  in 
further  details  of  the  various  investigations. 


12.1  HYPERVELOCITY  IMPACT 


The  selected  hypervelocity  impact  results  shown  here 
were  taken  from  a  numerical  parameter  study  [ 1 3 j  in  which 
material  properties  were  varied  in  order  to  assess  the  effects 
of  those  properties  on  cratering  and  mass  loss. 

Figures  12.1  and  12.2  show  the  predicted  projectile/ 
target  configuration  at  various  times  for  spheres  of  equal 
masses  of  ice  and  water  impacting  normally  into  an  aluminum 
target  at  a  velocity  of  3.05  x  10^  cm/sec. 

It  is  seen  from  Figure  12.3  that  there  is  no  significant 
difference  in  the  predicted  crater  depths  for  the  two  impacts 
involving  ice  and  water  projectiles.  The  fact  that  the 
maxima  in  the  crater  depth  curves  do  not  occur  at  the  latest 
times  is  a  result  of  crater  relaxation  which  is  typical  in 
soft  aluminum  targets.  Figure  12.4  is  a  comparison  of  the 
predicted  final  crater  shape  (solid  line)  for  the  ice- 
aluminum  impact  with  the  measured  profile  (dashed  lines) 
obtained  from  experimental  data  provided  by  Mr.  William 
Gray  at  Martin  Marietta  in  Orlando.  As  is  seen  in  Figure 
12.4,  agreement  between  theory  and  experiment  is  quite  good. 
The  shaded  region  in  the  vicinity  of  the  crater  lip  (Figure 
12.4)  shows  the  predicted  crater  ejecta  which  has  been 
removed  for  comparison  with  experiment  in  Figure  12.4. 

The  results  of  a  calculation  involving  a  spherical 
glass  projectile  impacting  6061-T6  aluminum  at  3.05  x  10^ 
cm/sec  was  also  compared  with  experimental  data  provided  by 
Martin  Marietta.  Figure  12.5  compares  the  calculated 
(solid  line)  and  experimentally  measured  (dashed  line) 
final  crater  profiles.  The  experimental  crater  was  not 
exactly  axisymmetric  so  profiles  taken  at  two  different  cross 
sections  are  shown  in  the  figure.  Both  the  predicted  crater 
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a)  Shaded  region  indicates  the  predicted  crater  ejecta. 
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b)  Final  predicted  crater  with  ejecta  removed. 

Figure  12 . 5- -Comparison  of  the  calculated  and  measured  craters 
for  a  0.953  cm  diameter  glass  sphere  impacting 
aluminum  at  3^05  x  10^  cm/sec. 
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depth  and  crater  diameter  are  in  good  agreement  with  ex¬ 
periment.  Again,  the  shaded  region  (Figure  12.5a)  has 
been  removed  in  Figure  12.5b  since  it  represents  failed 
material. 

Figure  12.6  shows  the  predicted  final  configurations 
for  the  calculation  of  a  glass  bead  into  ATJ-S  graphite  at 
6.1  x  10  cm/sec.  The  shaded  regions  represent  material 
failure.  In  addition  to  the  failure  predicted  in  the 
vicinity  of  the  crater  and  the  lip,  several  radial  bands  of 
failed  material  are  in  evidence.  This  type  of  failure  seems 
to  be  typical  of  ATJ-S  graphite  when  impacted  by  particles 
in  the  hypervelocity  regime. 
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12.2  PLUGGING  FAILURE 


When  a  blunt  projectile  impacts  a  relatively  thin 
target,  a  plug  of  target  material  is  sheared  out  ahead  of 
the  projectile.  Plugging  failure  is  modeled  in  HELP  by 
allowing  a  slip  surface  to  progress  through  the  target  in 
the  direction  of  the  maximum  shear  stress  when  the  material 
ahead  of  the  advancing  slipline  is  subjected  to  conditions 
such  that  the  plastic  work  exceeds  a  critical  value.  (See 
Chapter  VI.)  The  results  presented  in  this  subsection  were 
selected  from  several  recent  numerical  investigations  of 
plugging  failure  [14-17]. 

Figure  12.7  is  a  series  of  computer  output  plots  showing 
the  predicted  steel  projectile  and  aluminum  target  configura¬ 
tion  for  an  impact  calculation  at  various  times  after  impact. 
The  initial  impact  velocity  was  5.56  x  104  cm/sec.  Figure 
12.8  is  a  similar  plot  for  a  calculation  in  which  the  impact 
velocity  was  2.36  x  104  cm/sec.  Both  calculations  were 
carried  out  in  time  until  the  kinetic  energy  of  the  projec¬ 
tile  approached  its  asymptotic  value.  Several  differences 
in  the  computed  results  can  be  noted  by  comparing  Figures 

12.7  and  12.8.  The  plug  forms  much  sooner  in  the  higher 
velocity  impact  situation  (V0  =  5.56  x  104  cm/sec)  of 
Figure  12.7  than  in  the  lower  velocity  impact  situation 
(VQ  =  2.36  x  104  cm/sec)  of  Figure  12.8.  As  might  be  ex¬ 
pected,  Figure  12.7  indicates  that  the  projectile  deforms 
considerably  whereas  very  little  projectile  deformation  is 
evident  in  Figure  12.8.  Another  numerically  predicted 
velocity-related  trend  in  the  calculations  which  agrees 
with  observation  and  is  in  evidence  in  Figures  12.7  and 

12.8  involves  the  decrease  in  plug  taper  angle  as  the 
impact  velocity  decreases. 
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Figure  12.8--Projectile-target  configurations  at  variou 
times  for  the  calculation  of  a  steel  cylin 
lmpacting^an  aluminum  target  at  a  velocity 
2. 36  x  104  cm/sec. 
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cylinder  into  a  steel  target  at  a  velocity  of  8.81  x  104 
cm/sec.  The  effect  of  the  increased  target  density  and 
yield  strength  is  clearly  evident  by  noting  the  extensive 
projectile  deformation. 

A  plot  of  residual  projectile  velocity  versus  impact 
velocity  is  shown  in  Figure  12.10.  The  predicted  values 
from  these  calculations  (solid  circles)  are  compared  with 
experimental  data  (open  circles)  provided  by  Eglin  Air 
Force  Base  [16].  The  agreement  between  the  numerical 
predictions  and  experiment  is  quite  good. 
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The  plane  strain  option  in  HHLP  was  employed  to 
approximate  a  series  of  long  rod  oblique  impacts.  It  is 
important  to  keep  in  mind,  when  analyzing  the  results  from 
such  a  series  of  calculations,  that  the  real  three  dimensional 
case  is  not  being  solved.  For  example,  no  strain  is  allowed 
in  the  direction  perpendicular  to  the  plane  of  impact  (plane 
formed  by  initial  impact  velocity  vector  and  target  normal) . 

In  addition  the  mass  of  the  projectile  is  based  on  the 
volume  obtained  by  translating  the  central  cross  section  of 
the  projectile  one  unit  length  perpendicular  to  the  plane  of 
impact.  However,  many  things  can  be  learned  from  the  plane 
strain  approximation  to  oblique  impact.  In  a  recent  parameter 
study  [18]  it  was  of  interest  to  rank  several  different 
projectile  designs  according  to  the  extent  of  nose  tip  failure 
and  projectile  bending.  It  was  felt  that  these  rankings  could 
be  accomplished  by  comparing  results  from  a  series  of  2-D, 
plane-strain  calculations.  In  the  sample  calculation  presented 
here,  a  steel- jacketed,  tungsten  alloy  penetrator  is  impacted 
at  60°  obliquity  against  a  thin  steel  target  at  a  velocity  of 
1.53  x  10^  cm/sec.  The  pro j ecti le- target  configurations  at 
various  times  for  this  impact  solution  are  shown  in  Figure 


Figure  12 . 11 - -Two-dimensional  plane  strain  approximation  to 
an  oblique  impact. 


12.4  MULTIPLE  IMPACT 


The  HELP  code  was  employed  in  an  exploratory  investiga¬ 
tion  of  the  effect  of  projectile  geometry  on  back  surface 
spall  [19].  The  initial  configurations  of  the  two  impact 
situations  of  interest  are  shown  in  Figure  12.12.  The  solid 
projectile  is  a  steel  cylinder  with  a  hemispherical  nose. 

The  spaced  pellet  projectile  is  composed  of  four  steel 
pellets  having  a  total  mass  (3g) ,  equal  to  the  mass  of  the 
solid  projectile.  Both  projectiles  were  impacted  into  a 
1-cm-thick  steel  plate  at  a  velocity  of  9.15  x  10^  cm/sec. 

In  order  to  compute  the  impact  of  spaced  pellets  the  basic 
HELP  code  was  modified  so  that  voids  between  the  pellets 
could  be  removed  as  the  pellets  impacted  one  another. 

While  it  was  of  primary  interest  to  compare  the  two 
projectile  designs  in  terms  of  their  ability  to  create 
stress  conditions  suitable  for  back  surface  spall,  it  was 
also  of  interest  to  compare  the  penetrability  of  the  two 
projectiles.  Figures  12.13  and  12.14  are  computer  plots 
of  the  projectile- target  configurations  at  various  times 
after  impact  for  the  solid  and  spaced  pellet  projectiles, 
respectively.  The  crater  at  10  ysec  in  Figure  12.13  should 
be  compared  with  the  crater  at  12  ysec  in  Figure  12.14 
since  the  percent  reduction  in  kinetic  energy  for  the  two 
projectiles  is  equal  at  those  times.  This  time  lag  in 
projectile  kinetic  energy  reduction  is  a  result  of  the 
spacing  in  the  spaced-pellet  projectile.  It  can  be  concluded 
that  the  penetrability  of  the  two  projectile  designs  are 
equivalent  for  the  impact  situations  considered  here.  Situa¬ 
tions  could  exist,  however,  for  which  multiple  impact  might 
enhance  the  crater  size.  Such  enhancement  would  be  a  func¬ 
tion  of  relative  dimensions  and  pellet  spacing. 
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Figure  12 . 14- -Pro j ectile- target  configurations 
at  various  times  for  the  spaced 
pellet  projectile  impact. 


12.5  FRAGMENTATION  MUNITIONS 


This  subsection  presents  results  from  a  HELP  solution 
to  a  fragmentation  munitions  problem.  The  calculation  in¬ 
volved  a  hollow  steel  cylinder  filled  with  Octol.  Detona¬ 
tion  was  end- initiated  at  a  point  on  the  axis.  This  calcula¬ 
tion  was  taken  from  a  high  explosive  parameter  study  [20]. 
Figure  12.15  shows  the  initial  configuration  for  the 
calculation . 

The  high  explosive-casing  configurations  for  various 
times  after  detonation  initiation  are  shown  in  Figure  12.16. 
The  calculation  was  continued  in  time  until  failure  had  been 
predicted  throughout  the  entire  casing.  The  maximum  pressure 
in  the  detonated  products  by  that  time  had  been  reduced  to 
approximately  one  kilobar  and  the  casing  acceleration  was 
small. 

Figure  12.17  is  a  composite  plot  of  the  casing  at 
various  times  after  detonation  initiation.  The  shaded 
regions,  which  indicate  failed  material ,  show  that  failure 
occurs  first  at  the  outside  surface  of  the  casing  and  moves 
inward  until  the  entire  casing  has  failed.  In  the  calculation 
the  stresses  in  failed  cells  were  set  to  zero. 

Figure  12.18  gives  the  HELP  predictions  of  fragment 
velocity  and  projection  angle.  These  quantities  are  plotted 
versus  the  initial  relative  position  of  the  fragment  and  are 
compared  with  experimental  data  provided  by  the  BRL  [21]. 


12-22 


,,rlyr,  ’■ 


u . in,, , v,iu ,t  i v* v ' i/^ /. -».yv» 

i  ^  V *  ’■*.  ,_•*»  ■*  t  W  ta  .  *<  „  *\  *  -  *  -.  -•.».**..•-••*  *  **  -  "  *■  *  >  “  *  -’  *  ■*  "  * 


12.6  SHAPED  CHARGE  JET  FORMATION 

Recently  HELP  was  used  to  perform  several  numerical 
studies  involving  shaped  charge  jet  formation  [22-26].  The 
results  shown  here  are  selected  from  Reference  26  and 
involve  identical  copper  liners  (21°  half  angle)  and  Comp  B 
high  explosive,  but  different  casings.  Figure  12.19  shows 
the  initial  configurations  of  the  two  selected  calculations, 
and  Figure  12.20  and  12.21  show  the  predicted  configurations 
at  various  times  after  initiation  for  the  charges  confined 
in  aluminum  and  steel,  respectively.  Transmittive 
boundaries  at  the  bottom  of  the  calculational  grid  allow  the 
jet  material  to  leave  the  grid;  however,  quantities  of  in¬ 
terest  such  as  jet  mass  and  velocity  are  retained. 

Experimentally  determined  jet  tip  velocity  data  were 
made  available  by  the  BRL  [27]  for  comparison  with  the 
results  obtained  from  the  calculations  involving  aluminum 
confinement.  Figure  12.22  is  a  plot  of  calculated  jet 
velocities  versus  the  measured  jet  tip  velocities.  The 
open  circles  correspond  to  the  maximum  calculated  jet 
velocities.  These  calculated  maxima  did  not  occur  at  the 
jet  tip  because  of  the  inverse  velocity  gradient  in  that 
region,  but  they  do  fall  well  within  the  10%  error  limit 
shown  by  the  dashed  lines  in  Figure  12.22.  However  they 
are  consistently  slightly  lower  than  the  observed  values. 
This  is  probably  due  to  the  fact  that  numerically  calculated 
detonation  fronts  are  somewhat  broader  than  real  detonation 
fronts,  a  result  expected  unless  extremely  fine  resolution 
is  employed  in  the  calculations. 
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Figure  12 . 19- - Initial  configurations  of  two  calculations 

involving  copper  liners,  Comp  B  high  explosive 
and  light  (aluminum)  and  heavy  (steel)  casings 
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12.7  SHAPED  CHARGE  JET  PENETRATION 

The  HELP  code  has  also  recently  been  employed  in  an 
investigation  involving  shaped  charge  jet  penetration  of 
an  armor  plate  [28,  29].  The  calculations  were  carried  out 
in  time  until  the  jet  perforated  the  armor  and  the  region 
of  predicted  material  failure,  which  gives  rise  to  behind - 
the-armor  debris,  had  ceased  to  grow. 

The  initial  conditions  for  the  calculation  presented 
here  involve  a  copper  jet  having  a  bulb  at  the  tip,  impacting 
a  2.54  cm  thick  steel  plate  at  a  constant  velocity  of  7.6  x 
105  cm/sec.  The  bulb,  which  comprised  the  leading  3.81  cm 
of  the  jet,  was  described  by  a  straight  line  connecting  the 
1.27  cm  diameter  hemispherical  tip  to  the  constant  diameter 
portion  of  the  jet.  Figure  12.23  shows  the  predicted  jet  and 
armor  configurations  at  various  times  after  impact.  Figure 
12.24  shows  the  predicted  configuration  10.30  ysec  after 
initial  impact,  with  the  failed  cells  shaded.  The  total 
stress  components  were  set  to  zero  in  all  of  the  failed  cells. 

Some  of  the  results  from  the  calculation  can  be  compared 
with  available  experimental  data  obtained  from  the  BRL. 

Figure  12.25  is  a  comparison  of  the  numerically  predicted  debris 
cloud,  the  leading  edge  of  which  is  shown  at  various  times 
after  impact,  with  a  radiograph  of  actual  behind-the-armor 
debris.  The  numerical  prediction  agrees  well  with  the  ex¬ 
periment  . 
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Figure  12. 25- -Radiograph  of  behind- the -armor  debris  super¬ 
imposed  upon  the  predicted  shape  of  the  debris 
cloud  shown  at  various  times  after  impact. 
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12.8  WAVE  PROPAGATION  STUDY 

The  HELP  code  has  also  been  applied  to  the  study 
of  wave  propagation  in  various  media.  A  series  of  calcu¬ 
lations  was  performed  using  a  nitromethane  source  buried 
in  media  of  various  porosities  [30].  Stress-time  histories 
were  obtained  at  stations  located  at  several  ranges  from  the 
center  of  the  charge. and  compared  with  experimental  results 
[30].  Some  of  these  calculated  histories,  and  associated 
data  points,  are  presented  in  Figures  12.26  -  12.29.  The 
first  three  figures  were  obtained  at  shot  depth,  at  a  radius 
of  3  meters  from  the  center  of  the  source,  in  soils  of  in¬ 
creasing  porosity.  The  last  figure,  obtained  in  the  calcu¬ 
lation  of  the  soil  having  the  least  porosity,  is  the  history 

at  the  station  6  meters  from  the  center  of  the  source  at 
shot  depth. 

Figures  12.26  -  12.28  show  clearly  the  attenuation  of 
the  shock  strength  as  the  soil  porosity  increases.  Each  of 
the  figures  gives  excellent  agreement  with  the  data  as  regards 
shock  strength,  time  of  arrival,  and  the  shape  of  the  release 
curve.  Figure  12.29  shows  that  the  agreement  was  obtained 
even  at  stress  levels  on  the  order  of  1  kilobar. 

It  is  worth  noting  that  Figures  12.27  and  12.28  were 
calculated  as  part  of  a  design  study  prior  to  the  field 
experiment  from  which  the  data  were  obtained.  The  other 
two  figures  were  calculated  subsequent  to  the  running  of 
the  experiment. 


12-36 


ww 


Time  (msec) 


Figure  12. 26- -Predicted  stress  history  of  3  meter  station 
at  shot  depth  for  nitromethane  source  in 
1.3%  porosity  soil.  Data. is  from  a  test 
of  the  source  in  that  material. 
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Figure  12 . 27- -Predicted  stress  history  of  3  meter  station 
at  shot  depth  for  nitromethane  source  in  2% 
porosity  soil.  Data  is  from  a  test  of  the 


source  in  that  material. 
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Figure  12 . 28- -Predicted  stress  history  of  3  meter  station 
at  shot  depth  for  nitromethane  source  in 
4.2%  porosity  soil.  Data  is  from  a  test  of 
the  source  in  that  material. 
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Figure  12 . 29- -Predicted  stress  history  of  6  meter  station 
at  shot  depth  for  nitromethane  source  in 
1.3%  porosity  soil.  Data  is  from  two  different 
tests  of  the  source  in  that  material. 
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APPENDIX  A 


DIMENSIONING  THE  ARRAYS 


Dimensioning  the  arrays  of  the  HELP  code  to  enable  it 
to  handle  problems  of  arbitrary  size  (limited  only  by  the 
core  memory  of  the  user's  computer)  is  a  relatively  straight¬ 
forward  procedure.  It  is  necessary  only  to  specify  16 
primary  parameters  and  22  secondary  parameters  in  order  to 
uniquely  determine  all  array  dimensions  and  equivalences  in 
the  HELP  code. 

The  primary  parameters  are  listed  in  Table  Al,  along 
with  a  description  of  each  and,  when  it  exists,  the  name  of 
the  corresponding  HELP  variable.  Because  of  equivalencing 
between  various  arrays,  there  are  constraints  placed  on  the 
values  assigned  to  some  of  these  parameters.  These 
constraints  are  listed  in  Table  A2. 

After  the  user  has  determined  the  values  of  the  16 
primary  parameters,  it  is  necessary  to  determine  another  22 
secondary  parameters,  which  are  functions  of  the  primary 
parameters.  These  secondary  parameters  are  defined  in 
Table  A3. 

Having  defined  these  parameters,  all  that  remains  is 
for  the  user  to  incorporate  their  values  into  the  common, 
dimension,  and  equivalence  statements  that  use  parameter 
rather  than  integer  dimensions.  Most  of  these  statements 
are  in  the  "included"  element,  and  the  remaining  ones  are  in 
subroutine  SPHASE.  Any  arrays  which  are  dimensioned  by 
integer  values  (not  parameters)  are  fixed  and  should  not 
be  redimensioned. 


TABLE  A . 1 


PRIMARY  PARAMETERS  POR  DIMENSIONING  THE  HELP  CODE 


Parameter 

Definition* 

MDX 

Maximum  number  of  interface  cells  which  can 

occur  on  any  one  cycle.  (NMXCLS) 

NMX 

Number  of  material  packages  in  the  problem, 
excluding  the  void  package.  (NMAT) 

NTPX 

Maximum  number  of  material  tracer  particles  for 
any  one  package  in  the  problem.  (NTPMX) 

IDX 

Number  of  columns  in  the  grid.  (IMAX) 

See  Table  A2. 

JDX 

Number  of  rows  in  the  grid.  .  (JMAX)  See  Table  A2. 

I JDX 

If  problem  is  to  be  rezoned,  IJDX  is  the  maximum 
value  of  IDX  and  JDX.  Otherwise,  it  can  be  set 

to  1. 

KSTR 

If  problem  will  compute  strength  effects 
CCYCPH3>0)  ,  set  KSTR  to  IDX^JDX  +  1.  Otherwise, 

set  to  1. 

IDXX 

If  problem  will  compute  strength  effects,  set 

IDXX  to  IDX  +  1.  Otherwise  set  to  1.  This 

parameter  is  used  only  to  dimension  local 

variables  in  SPHASE. 

The  analagous  HELP  variable,  if  any,  is  in  parentheses 


A-  2 


Parameter 


Def ini tion 


KDET  If  problem  will  detonate  an  explosive,  set  KDET 

to  IDX  •  JDX  +  1.  Otherwise,  set  to  1, 

KPX  Dimension  of  pressure  array.  See  Table  A. 2  for 

determining  the  value  of  this  parameter. 

NSLDE  If  the  problem  will  have  a  slip  line  (NOSLIP=0) , 

set  NSLDE  to  the  maximum  number  of  mixed  cells 
which  will  contain  the  slip  line.  Otherwise, 
set  it  to  1.  (NSLD) 

JDX2  Dimension  of  UL  and  PL  arrays.  See  Table  A. 2  for 

determining  the  value  of  this  parameter. 

IPLGR  If  problem  is  a  plugging  calculation,  set  IPLGR 

to  the  value  of  IPLGRT,  defined  in  Section  7.2.1. 
Otherwise  set  to  1.  (IPLGRT) 

IPLGZ  If  problem  is  a  plugging  calculation,  set  IPLGZ 

to  the  number  of  rows  in  the  target.  Otherwise 
set  to  1.  (IPLGTP-IPLGBT  +  1.)  See  definition 
of  IPLGTP  and  IPLGBT  in  Section  7.2.1. 

KDX4  Maximum  number  of  passive  tracer  particles  the 

problem  will  generate.  If  the  problem  is  a 
plugging  calculation,  set  KDX4  to  4  • IPLGR  •  IPLGZ . 
If  passive  tracers  are  not  generated  (NTCC  =  0) , 
set  KDX4  to  1.  (NTCC) 

K4PL  Dimension  of  PLWP  array.  If  problem  is  a  plugging 

calculation,  set  K4PL  to  4  •  IPLGR  •  IPLGZ . 

Otherwise  set  to  1. 


TABLE  A. 2 


CONSTRAINTS  ON  PARAMETER  VALUES 


Parameter 


Applicable  Constraint 


IDX 

JDX 

KPX 


Must  be  at  least  3. 

Must  be  at  least  3. 

Set  to  maximum  value  of  the  following 


IDX  •  JDX  +  1 
2  •  (NMX  +  1)  •  MDX 
NMX  •  (5  •  JDX  +12) 

10  •  NMX  +  (4  +  3  •  NMX)  • 


JDX2 


Set  to  maximum  value  of  the  following 
IDX 


I JDX  +  2 


TABLE  A. 3 


APPENDIX  B 


ABBREVIATED  HELP  INPUT  INSTRUCTIONS 

After  becoming  familiar  with  the  various  input  options 
and  specifications  discussed  in  Section  7.2,  the  user  can 
use  the  following  forms  to  quickly  describe  the  input  for 
a  given  HELP  calculation.  The  default  values  of  most 
Z-block  variables  are  given  in  parentheses  after  the 
definition.  If  no  value  is  specified,  the  default  value 
is  zero. 


B  - 1 


The  irput  variables  are  described  in  more  detail  in  Chapter  VII.  Instructions 
restarting  a  calculation  are  given  in  Section  7.3. 
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APPENDIX  C 


SEGMENTING  THE  HELP  CODE 

The  user  of  the  HELP  code  has  several  alternatives  in 
configuring  the  code  on  the  computer,  ranging  from  having 
all  the  subroutines  loaded  into  core  at  once,  which  decreases 
execution  time  but  requires  more  core  storage,  to  varying 
degrees  of  segmentation,  which  decrease  the  core  storage 
requirements  but  increase  execution  time  because  of  the 
extra  time  required  to  transfer  the  various  segments  in  and 
out  of  core.  The  segmentation  discussed  here  and  represented 
schematically  in  Figure  C.l  attempts  to  minimize  the  storage 
requirements,  yet  hold  execution  time  down  by  minimizing  the 
number  of  times  per  cycle  the  computer  must  roll  segments  in 
and  out  of  core. 

In  the  figure,  the  subroutines  in  the  first  column, 
which  is  labeled  LEVEL  1,  are  loaded  into  core  and  remain  in 
core  for  the  duration  of  the  calculation.  The  blocks  of  sub¬ 
routines,  or  segments,  under  LEVEL  2  and  LEVEL  3  are  then 
rolled  in  and  out  of  core  as  needed  by  the  various  stages 
of  the  computational  cycle.  The  INPUT -GENERATOR  segment  is 
executed  only  once  per  run,  either  in  generating  or  re¬ 
starting  a  calculation.  If  the  user  is  not  performing  a 
plugging  calculation,  the  REZONE-PLUG  block  is  executed  only 
when  the  grid  is  rezoned,  a  relatively  infrequent  occurence. 

Thus  in  most  of  the  cycles  only  four  LEVEL  2  segments  (the 
CDT-EDIT  segment,  the  HPHASE- SPHASE  segment,  the  INFACE 
segment,  and  the  TPHASE  segment)  and  two  LEVEL  3  segments 
are  rolled  in  and  out  of  core. 

This  configuration  of  the  HELP  code,  and  the  associated 
systems  routines,  requires  approximately  11,000  words  of  core 
storage  for  the  instruction  bank  when  executed  on  a  UNIVAC  1108. 
The  amount  of  storage  required  for  the  data  bank  is,  of  course 
highly  dependent  on  the  values  of  the  parameters  discussed  in 
Appendix  A. 
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